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Abstract 

Future  applications  in  environmental  monitoring,  delivery  of  services  and  transportation  of  goods  moti¬ 
vate  the  study  of  deployment  and  partitioning  tasks  for  groups  of  autonomous  mobile  agents.  These  tasks 
are  achieved  by  recent  coverage  algorithms,  based  upon  the  classic  methods  by  Lloyd.  These  algorithms  how¬ 
ever  rely  upon  critical  requirements  on  the  communication  network:  information  is  exchanged  synchronously 
among  all  agents  and  long-range  communication  is  sometimes  required.  This  work  proposes  novel  coverage 
algorithms  that  require  only  gossip  communication ,  i.e.,  asynchronous,  pairwise,  and  possibly  unreliable 
communication.  Which  robot  pair  communicates  at  any  given  time  may  be  selected  deterministically  or 
randomly.  A  key  innovative  idea  is  describing  coverage  algorithms  for  robot  deployment  and  environment 
partitioning  as  dynamical  systems  on  a  space  of  partitions.  In  other  words,  we  study  the  evolution  of  the 
regions  assigned  to  each  agent  rather  than  the  evolution  of  the  agents’  positions.  The  proposed  gossip 
algorithms  are  shown  to  converge  to  centroidal  Voronoi  partitions  under  mild  technical  conditions. 

Our  treatment  features  a  broad  variety  of  results  in  topology,  analysis  and  geometry.  First,  we  establish 
the  compactness  of  a  suitable  space  of  partitions  with  respect  to  the  symmetric  difference  metric.  Second, 
with  respect  to  this  metric,  we  establish  the  continuity  of  various  geometric  maps,  including  the  Voronoi 
diagram  as  a  function  of  its  generators,  the  location  of  a  centroid  as  a  function  of  a  set,  and  the  widely-known 
multicenter  function  studied  in  facility  location  problems.  Third,  we  prove  two  convergence  theorems  for 
dynamical  systems  on  metric  spaces  described  by  deterministic  and  stochastic  switches. 


1  Introduction 

In  the  not  too  distant  future,  networks  of  coordinated  autonomous  robots  will  perform  a  broad  range  of  en¬ 
vironmental  monitoring  and  logistic  tasks.  Robotic  camera  networks  will  monitor  airports  and  other  public 
infrastructures.  Teams  of  vehicles  will  perform  surveillance,  exploration  and  search  and  rescue  operations. 
Groups  of  robots  will  enable  novel  logistic  capacities  in  the  transportation  of  goods  and  the  delivery  of  services 
and  resources  to  customers.  New  applications  will  be  enabled  by  the  ongoing  decreases  in  size  and  cost  and  the 
increases  in  performance  of  sensors,  actuators,  communication  devices  and  computing  elements. 

In  these  future  applications,  load  balancing  algorithms  will  dictate  how  the  workload  is  shared  amongst  and 
assigned  to  the  individual  robots.  In  other  words,  robotic  resources  will  be  assigned  and  deployed  to  competing 
requests  in  such  a  way  as  to  optimize  some  performance  metric.  Remarkably,  load  balancing  problems  in  robotic 
networks  are  often  equivalent  to  robotic  deployment  and  environment  partitioning  problems.  For  example,  in 
surveillance  applications,  optimal  sensor  coverage  is  often  achieved  by  partitioning  the  environment  and  assigning 
individual  robotic  sensors  to  individual  regions  of  responsibility.  Similarly,  in  the  transportation  of  goods  or 
delivery  of  services,  minimizing  the  customer  wait-time  is  equivalent  to  a  multi- vehicle  routing  problem  and,  in 
turn,  to  computing  optimal  depot  positions  and  regions  of  responsibility. 
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Motivated  by  these  scenarios,  this  paper  considers  the  two  following  interrelated  problems.  The  deployment 
problem  for  a  robotic  network  amounts  to  the  design  of  coordination  algorithms  that  lead  the  robots  to  be  opti¬ 
mally  placed  in  an  environment  of  interest.  Deployment  performance  is  characterized  by  an  appropriate  network 
utility  function  that  measures  the  deployment  quality  of  a  given  configuration.  The  partitioning  problem  is  the 
design  of  coordination  algorithms  that  lead  the  robots  to  optimally  partition  the  environment  into  subregions 
of  interest;  even  here  the  objective  is  usually  achieved  through  the  design  of  appropriate  utility  functions. 

Literature  review 

A  wide  range  of  literature  is  relevant  to  our  investigation.  In  next  four  paragraphs  we  review  previous  work  on 
centroidal  Voronoi  partitions,  distributed  robotic  deployment,  animal  territory  partitioning,  and  mathematical 
modeling  tools. 

The  “centering  and  partitioning”  algorithm  originally  proposed  by  Lloyd  [25]  and  elegantly  reviewed  in  the 
survey  [T4]  is  a  classic  approach  to  facility  location  and  environment  partitioning  problems.  The  Lloyd  algorithm 
computes  centroidal  Voronoi  partitions  as  optimal  configurations  of  an  important  class  of  objective  functions, 
called  multicenter  functions.  Besides  their  intended  application  to  quantization  theory  he  centroidal  Voronoi 
partitions  have  widespread  applications  in  numerous  disciplines,  including  statistical  pattern  recognition  |22j. 
mathematical  imaging  HE  geometric  optimization  [3]  and  spatial  resource  allocation  [12]  ■  Recent  mathematical 
interest  has  focused  on  convergence  analysis  [13] .  bifurcation  analysis  of  low  dimensional  problems  [38] .  and 
anisotropic  partitions  HU,  among  other  topics. 

Distributed  and  robotic  versions  of  the  Lloyd  algorithm  have  been  recently  developed  in  the  multiagent 
literature;  see  the  survey  [28]  and  the  text  [9]  Chapter  5  and  literature  notes  in  Section  5.4],  We  briefly 
review  this  growing  literature  in  what  follows.  Generalized  centroidal  Voronoi  partitions  are  shown  in  [171  to 
be  asymptotically  optimal  for  estimation  of  stochastic  spatial  fields  by  sensor  networks.  In  m  convergence 
results  are  obtained  for  simple  basic  interactions  among  mobile  agents  such  as  “move  away  from  the  closest 
other  agent”  or  “move  toward  the  furthest  vertex  of  your  region  of  responsibility.”  Convergence  to  centroidal 
Voronoi  partitions  is  established  in  [4]  for  a  class  of  communication-less  sensor-based  algorithms  (related  to  the 
classic  clustering  work  by  MacQueen  [27] ).  In  [36]  adaptive  coverage  controls  are  proposed  for  environments 
described  by  unknown  density  functions.  Dynamic  environments  and  corresponding  dynamic  coverage  problems 
are  treated  in  [21! .  In  [31]  partitioning  policies  are  shown  to  achieve  optimal  load  balancing  in  vehicle  routing 
problems,  i.e. ,  problems  in  which  a  robotic  network  provides  service  to  customers  that  arrive  in  real  time  in  the 
environment. 

Territory  partitioning  via  competitive  behaviors  is  a  classic  subject  of  study  in  behavioral  ecology;  see  [I] 
for  a  comprehensive  survey.  For  example,  it  is  known  [2]  that  the  foraging  behavior  of  conflicting  colonies 
of  red  harvester  ants  ( Pogonomyrmex  barbatus )  results  in  non-overlapping  dominance  regions  that  resemble 
Voronoi  partitions.  Non-overlapping  dominance  regions  akin  to  centroidal  Voronoi  partitions  are  well  docu¬ 
mented  in  Ensma  for  the  mouthbreeder  fish  ( Tilapia  mossambica).  Territoriality  behavior  and  competition 
among  prides  of  African  lions  ( Panthera  leo)  are  discussed  in  [29] .  Overall,  numerous  animal  species  achieve  ter¬ 
ritory  partitioning  without  a  central  coordinating  entity  and  without  synchronized  communication,  but  rather 
relying  upon  asynchronous  accidental  interactions  and  stigmergy.  To  the  best  of  our  knowledge  about  biological 
and  engineering  multiagent  systems,  asynchronous  territory  partitioning  has  been  barely  studied,  see  [37]  for 
introductory  ideas  about  animal  behavior,  and  mathematical  models  and  analysis  are  lacking. 

To  finalize  the  literature  review,  here  is  a  synopsis  of  the  mathematical  ideas  from  diverse  disciplines  that 
we  bring  to  bear  on  deployment  and  partitioning  problems.  First,  we  adopt  the  so-called  gossip  communication 
model ,  in  which  only  peer-to-peer  asynchronous  communication  links  are  required.  This  communication  model 
is  widely  studied  in  the  wireless  communication  literature;  example  references  include  [23]  [7].  Moreover,  we 
consider  control  systems  on  a  non-Euclidean  state  space;  the  interest  for  non-Euclidean  spaces  has  a  rich  history 
in  nonlinear  control  theory,  dynamical  systems  and  robotics,  including  the  early  work  [8]  and  a  recent  application 
to  multiagent  systems  |35j.  Finally,  we  adopt  various  tools  from  topology  and  from  the  study  of  hyperspaces  of 
sets  [30]. 
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Statement  of  contributions 

This  paper  uncovers  novel  mathematical  principles  and  tools  of  relevance  in  coordination  problems  and  multi¬ 
agent  systems. 

We  tackle  partitioning  and  coverage  control  algorithms  in  innovative  ways.  First,  we  design  algorithms  that 
require  only  gossip  communication,  i.e.,  asynchronous,  pairwise,  and  possibly  unreliable  communication.  Gossip 
communication  is  a  simple,  robust  and  effective  protocol  for  noisy  and  uncertain  wireless  environments.  Gossip 
communication  may  be  implemented  in  wandering  robots  with  short-range  unreliable  communication  (see  the 
illustrative  motion  coordination  strategy  in  Appendix  [Aj.  Second,  we  propose  a  paradigm  shift  in  coverage 
control  and  multicenter  optimization.  Classically  da  El  EH  H2  EU  H2  EHj,  the  state  space  for  the  coverage 
algorithms  are  the  agents’  positions,  i.e.,  as  a  function  of  the  agents’  positions  the  environment  is  divided  into 
regions  and  regions  are  assigned  to  each  agent.  Instead,  in  our  approach,  the  agents’  positions  are  no  longer  a 
concern:  the  state  space  is  the  space  of  partitions  of  the  environment  and  the  algorithm  dictates  how  to  update 
the  regions. 

Within  the  innovative  context  of  gossip  communication  and  partition-based  mechanisms,  we  devise  a  novel 
algorithm  for  multicenter  and  coverage  optimization.  Our  gossip  coverage  algorithm  is  a  peer-to-peer  version 
of  Lloyd  algorithm  and  aims  to  compute  centroidal  Voronoi  partitions.  Which  robot  pair  communicates  at 
any  given  time  is  the  outcome  of  either  a  deterministic  or  a  stochastic  process.  We  also  propose  a  modified 
version  that  restricts  communication  exchanges  to  adjacent  regions  and  that  has  suitable  continuity  properties. 
Simulations  illustrate  that  our  algorithms  successfully  compute  centroidal  Voronoi  partitions. 

To  formally  establish  the  convergence  properties  of  our  proposed  gossip  coverage  algorithms,  we  perform  a 
detailed  mathematical  analysis  composed  of  three  steps.  First,  we  develop  suitable  versions  of  the  Krasovskii- 
LaSalle  invariance  principle  for  dynamical  systems  on  metric  spaces  described  by  deterministic  and  stochastic 
switches.  Convergence  to  a  set  of  fixed  points  is  achieved  under  a  uniform  deterministic  or  stochastic  persistency 
condition.  Second,  we  establish  the  continuity  of  various  geometric  maps,  including  (1)  the  Voronoi  diagram  as 
a  function  of  its  generators,  (2)  the  location  of  the  generalized  centroid  as  a  function  of  a  set,  (3)  the  widely- 
known  multicenter  function  studied  in  facility  location,  and  (4)  the  gossip  coverage  algorithms.  These  continuity 
properties  are  established  with  respect  to  the  symmetric  distance  metric  in  the  space  of  partitions.  Third  and 
final,  we  study  the  topology  of  the  space  of  partitions.  With  respect  to  the  symmetric  difference  metric  we 
prove  the  compactness  of  a  relevant  subset  of  partitions.  Specifically,  we  focus  on  partitions  whose  component 
regions  are  the  union  of  a  bounded  number  of  convex  sets. 

In  summary,  relying  upon  our  extensions  of  the  invariance  principle,  the  compactness  of  a  subset  of  the  set 
of  partitions,  and  the  continuity  of  the  various  relevant  maps,  we  establish  the  convergence  properties  of  the 
proposed  gossip  algorithms.  In  short,  the  algorithms  converge  to  the  set  of  centroidal  Voronoi  partitions  under 
mild  technical  assumptions  and  under  the  assumption  that  the  gossip  communication  exchanges  satisfy  either 
a  deterministic  or  a  stochastic  persistency  condition. 

Organization  and  notations 

The  paper  is  structured  as  follows.  In  Section [2]  we  review  multicenter  optimization  and  coverage  control  ideas. 
In  Section[3]we  state  our  asynchronous  territory  partitioning  problem,  provide  a  solution  via  the  gossip  coverage 
algorithm,  state  the  convergence  properties  of  the  algorithm  and  report  some  simulation  results.  The  following 
Sections  U  [5]  and  [6]  develop  the  mathematical  machinery  required  to  prove  the  convergence  results.  Section  [4] 
contains  the  convergence  theorems  extending  the  Krasovskii-LaSalle  invariance  principle.  Section  [5]  contains  a 
discussion  about  the  compactness  properties  of  the  space  of  partitions.  Sect  ion  [6]  states  the  continuity  properties 
of  the  relevant  maps  and  functions  and  contains  the  proof  of  the  main  convergence  results.  Concluding  remarks 
are  given  in  Section  [7] 

We  let  R>o  and  R>o  denote  the  set  of  positive  and  non-negative  real  numbers,  respectively,  and  Z>0  denote 
the  set  of  non-negative  integer  numbers.  Given  A  C  Rd,  we  let  int(A),  A,  dA  and  diam(A)  denote  its  interior, 
its  closure,  its  boundary  and  its  diameter,  respectively.  Given  two  sets  X  and  V,  a  set-valued  map  T  :  X  =4  Y 
associates  to  an  element  of  A  a  subset  of  Y . 
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2  A  review  of  multicenter  optimization  and  distributed  coverage 
control 

In  this  section  we  review  a  variety  of  known  results  in  geometric  optimization  and  in  robotic  coordination.  In 
Subsection  I2J1  we  review  the  notion  of  environment  partitions  and  we  introduce  the  multicenter  function  as 
a  way  to  define  optimal  environment  partitions  and  optimal  robot  or  sensor  positions  in  the  environment.  In 
Subsection  W2\  we  review  some  distributed  control  algorithms  for  agent  motion  coordination  and  environment 
partitioning  based  on  the  classic  Lloyd  algorithm. 

2.1  Partitions,  centroids  and  multicenter  optimization 

We  let  Q  denote  an  environment  of  interest  to  be  apportioned.  We  assume  Q  is  a  compact  convex  subset  of 
with  non-empty  interior.  Partitions  of  Q  are  defined  as  follows. 

Definition  2.1  (Partition)  An  TV-partition  of  Q,  denoted  by  v  =  is  an  ordered  collection  of  N  subsets 

of  Q  with  the  following  properties: 

(i)  Vi  Q, 

(ii)  int(i’j)  fl  int(uj)  is  empty  for  all  i,j  G  {1, . . . ,  N}  with  i  ^  j ;  and 
(in)  each  set  Vi,  i  G  (1, . . .  ,IV},  is  closed  and  has  non-empty  interior. 

The  set  of  N -partitions  of  Q  is  denoted  by  Vat. 

Let  p  =  (pi, . . .  ,pn)  G  Qn  denote  the  position  of  N  agents  in  the  environment  Q.  Given  a  group  of  N  agents 
and  an  TV-partition,  each  agent  is  naturally  in  one-to-one  correspondence  with  a  component  of  the  partition; 
specifically  we  refer  to  Vi  as  the  dominance  region  of  agent  i  G  {1,...,1V}. 

On  Q,  we  define  a  density  function  to  be  a  bounded  measurable  positive  function  (f>  :  Q  — >  R>o  and  a 
performance  function  to  be  a  locally  Lipschitz,  monotone  increasing  and  convex  function  /  :  R>o  — >  R>o-  With 
these  notions,  we  define  the  multicenter  function  Hmuiticenter  :  Vat  x  Qn  — >  M>o  by 

N  f 

7"fmulticenter(G  p)  —  f(\\Pi-<l\\)cl>(q)dq.  (1) 

i=i  Jv * 

This  function  is  well-defined  because  closed  sets  are  measurable.  We  aim  to  minimize  Wmuiticenter  with  respect 
to  both  the  partition  v  and  the  locations  p. 

Remarks  2.2  (Locational  optimization)  As  discussed  in  the  introduction  and  in  the  survey  \14V,  the  mul¬ 
ticenter  function  7tmuitiCenter  has  numerous  interpretations.  Here  we  review  two  applications  entailing  robotic 
networks.  First,  in  a  vehicle  routing  and  service  delivery  example  m,  given  vehicles  at  locations  pt,  assume 
that  f{\\pi  —  ?||)  is  the  cost  incurred  by  agent  i  to  travel  to  service  an  event  taking  place  at  point  q.  Events 
take  place  inside  Q  with  likelihood  (j).  Accordingly,  7fmuitiCenter  quantifies  the  expected  wait-time  between  event 
arrivals  and  agents  servicing  them. 

Second,  in  an  environmental  monitoring  application  \28'f.  assume  the  robots  aim  to  detect  acoustic  signals 
that  originate  and  propagate  isotropically  in  the  environment.  Because  of  noise  and  loss  of  resolution,  the  ability 
to  detect  a  sound  source  originating  at  a  point  q  from  a  sensor  at  position  pi  is  proportional  to  the  signal-to-noise 
ratio  (which  degrades  with  ||g  —  Pi\\).  If  the  performance  function  f  equals  minus  the  signal-to-noise  ratio,  then 
Tfmuiticenter  quantifies  the  expected  signal-to-noise  ratio  and  detection  capacity  for  acoustic  signals  generated  at 
random  locations.  □ 

Among  all  possible  ways  of  partitioning  a  subset  of  M.d,  one  is  worth  of  special  attention.  Define  the  set  of 
partly  coincident  locations  Sn  =  {p  G  QN  \  Pi  =  Pj  for  some  i,j  G  { 1 , . . . ,  N},  i  j}.  Given  p  G  QN  \  Sn, 
the  Voronoi  partition  of  Q  generated  by  p,  denoted  by  V (p) ,  is  the  ordered  collection  of  the  Voronoi  regions 

(vi(p))iLv  defined  by 

Vfip)  =  {qeQ  \  ||(? -Pi||  <  Ik  “Pill  for  ad  j  ^  *}■  (2) 


4 


In  other  words,  the  Voronoi  partition  is  a  map  V  :  ( QN  \  Sn)  — >  Vn-  The  regions  Vi(p),  i  £  {1, . . . ,  IV},  are 
convex  and,  if  Q  is  a  polytope,  they  are  polytopes.  Now,  given  two  distinct  points  q\  and  q2  in  R  ,  define  the 
(91;  92) -bisector  half-space  by 


^bisector (91 5  92)  =  {9  e  Rd  I  II 9  -  91 1|  <  ||9  -  92II}-  (3) 

In  other  words,  the  set  Fbisector(9i;  92)  is  the  closed  half-space  containing  qi  whose  boundary  is  the  hyperplane 
bisecting  the  segment  from  91  to  92-  Note  that  bisector  subspaces  satisfy  -ffbisector(9i;  92)  7^  Fbisectorfe;  9i)  and 
that  Voronoi  partition  of  Q  satisfies  Vi(pi, . . .  ,pn)  =  Q  H  (  flj^i  -ffbisector {Pi',Pj))  ■ 

Each  region  equipped  with  a  density  function  possesses  a  point  with  a  special  relationship  with  the  multi¬ 
center  function.  Define  the  scalar  1  -center  function  Hi  by 

Wi(p;  A)  =  J  f(\\p-  q\\)c/>(q)dq,  (4) 

where  p  is  any  point  in  Q  and  A  is  a  compact  subset  of  Q.  Under  the  stated  assumptions  on  the  performance 
function  /,  the  function  p  1— >  Hi(p;  A)  is  strictly  convex  in  p,  for  any  set  A  with  positive  measure  (we  postpone 
the  proof  to  Lemma  16.111.  Therefore,  the  function  p  1— >  Hi(p\A)  has  a  unique  minimum  in  the  compact  and 
convex  set  Q.  We  define  the  generalized  centroid  of  a  compact  set  A  C  Q  with  positive  measure  by 

Cd(.A)  =  argmin{7U(p;  A)  \  p  £  Q}.  (5) 

In  what  follows,  it  is  convenient  to  drop  the  word  “generalized,”  and  to  denote  by  Cd(i>)  =  (Cd(iq), . . . ,  Cd(ujv))  € 
Qn  the  vector  of  regions  centroids  corresponding  to  a  partition  v  £  Vjv. 

Remark  2.3  (Quadratic  and  linear  performance  functions)  If  the  performance  function  is  f(x)  =  x2 , 
then  the  global  minimum  of  Hi  is  the  centroid  ( also  called  the  center  of  mass)  of  A,  defined  by 

Cd {A)  =  (^  J  (/)(q)dq^  J  q(l>(q)dq. 

If  the  performance  function  is  fix)  =  x,  then  the  global  minimum  of  Hi  is  the  median  (also  called  the  Fermat- 
Weber  center)  of  A.  See  fH  Chapter  2]  for  more  details.  □ 

Voronoi  partitions  and  centroids  have  useful  optimality  properties  stated  in  the  following  proposition  and 
illustrated  in  Fig.  |TJ 

Proposition  2.4  (Properties  of  7^muiticenter)  For  any  partition  v  £  Vn  and  any  point  set  p  £  QN  \  Sn, 


!"fmulticenter(  V(p)  ,  p)  V  I^multicenter(^?  p)  ,  (6) 

multicenter  (V  ^d(u))  multicenter (u  P)  ■  (7) 

Furthermore,  inequality  ([6|  is  strict  if  any  entry  of  V ( p )  differs  from  the  corresponding  entry  of  v  by  a  set  with 
positive  measure,  and  inequality  0  is  strict  if  Cd(u)  differs  from  p. 

The  statements  in  Proposition  12.41  originate  in  the  early  work  by  S.  P.  Lloyd  [25|:  modern  treatments  are 
given  in  m  and  10  Propositions  2.14  and  2.15].  Proposition  12.41  implies  the  following  necessary  condition:  if 
a  pair  (v,p)  with  p  £  Sn  minimizes  H multicenter;  then  p  =  Cd(u)  and  v  =  V(p)  up  to  a  set  of  measure  zero. 
Accordingly,  the  partitions  that  minimize  ?dmuiticenter  have  the  following  property. 

Definition  2.5  The  partition  v  £  Vn  is  centroidal  Voronoi  if  it  has  distinct  centroids,  that  is,  Cd(uj)  Cd(uj) 
for  all  j  i,  and  v  =  V(Cd(u)). 
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Figure  1:  Illustration  of  Proposition^. 41  The  left  figure  shows  a  sample  2-partition  v  and  point  set  p  in  a  uniform 
square  environment.  The  value  of  the  cost  function  Wmuiticenter  at  (v,p)  is  diminished  by  either  replacing  v  with 
the  Voronoi  partition  generated  by  p  (see  central  figure),  or  replacing  p  with  the  centroids  of  v  (see  right  figure). 


Figure  2:  The  Voronoi  partition  and  the  corresponding  Delaunay  graph 


2.2  From  Lloyd  algorithm  to  distributed  coverage  control 

Here  we  consider  a  group  of  robotic  agents  with  motion,  communication  and  computation  capacities  and  we 
review  a  coverage  control  algorithm  that  determines  the  motion  of  each  robot  in  a  group  and  the  associated 
partition  in  such  a  way  as  to  minimize  7tmuiticenter-  In  what  follows,  we  restrict  our  attention  to  d  =  2,  that  is, 
we  assume  Qcl2. 

To  explain  in  what  sense  our  algorithms  are  distributed,  we  introduce  a  useful  graph.  The  Delaunay  graph  PH 
[9]  associated  to  the  distinct  positions  p  £  QN  \  Sn  is  the  undirected  graph  with  node  set  {Pi}f=  i  and  with  the 
following  edges:  ( Pi,Pj )  is  an  edge  if  and  only  if  V) (p)  fl  Vj (p)  is  non-empty.  In  other  words,  two  agents  are 
neighbors  if  and  only  if  their  Voronoi  regions  intersect,  see  Fig.  EJ 

The  coverage  algorithm  we  consider  is  a  distributed  version  of  the  classic  Lloyd  algorithm  |14|  based  on 
“centering  and  partitioning”  for  the  computation  of  centroidal  Voronoi  partitions.  The  algorithm  is  distributed 
in  the  sense  that  each  robot  determines  its  region  of  responsibility  and  its  motion  plan  based  upon  communication 
with  only  some  neighbors.  Specifically,  communications  among  the  robots  takes  place  along  the  edges  of  the 
Delaunay  graph.  The  distributed  coverage  algorithm  is  described  as  follows.  At  each  discrete  time  instant 
t  £  Z>o,  each  agent  i  performs  the  following  tasks:  (i)  it  transmits  its  position  and  receives  the  positions  of 
its  neighbors  in  the  Delaunay  graph;  (ii)  it  computes  its  Voronoi  region  with  the  information  received;  (iii)  it 
moves  to  the  centroid  of  its  Voronoi  region.  In  mathematical  terms,  for  t  £  Z>o, 

p(t  +  l)  =  Cd(V(p(t))).  (8) 

A  variation  of  the  function  7YmuitiCenter  is  useful  to  analyze  this  algorithm.  We  define  the  positions-based 
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multicenter  function  Tfvt 


oronoi  • 


,i  :Qn\Sn^  R>0  by 


N 

Voronoi  (p)  =  7tmulticenter(V(p) 7  V)  ''f  '  /  /(||<7  Pi  II 

i=l  ^(p) 


oronoi 


(9) 


Because  of  the  compactness  of  Q,  a  continuity  property,  and  the  monotonicity  properties  in  Proposition^. 41  one 
can  show  El  Theorem  5.5]  that  Ttvoronoi  is  monotonically  non-increasing  along  the  solutions  of  ©  and  that  all 
solutions  of  ©  converge  asymptotically  to  the  set  of  configurations  that  generate  centroidal  Voronoi  partitions. 
Additional  considerations  about  convergence  are  given  in  m- 

3  Gossip  coverage  control  as  a  dynamical  system  on  the  space  of 
partitions 

In  this  section  we  present  the  problem  of  interest,  our  novel  gossip  coverage  algorithm  and  its  convergence 
properties  in  Subsections  13.11 13.21  and  13.31  respectively.  In  order  to  reduce  the  communication  requirements  of 
our  algorithm,  we  propose  an  adjacency-based  and  continuous  algorithm  in  Subsection  13.41  Finally,  we  report 
some  simulation  results  in  Subsection  1331 

3.1  Problem  statement 

The  distributed  coverage  law,  based  upon  the  Lloyd  algorithm  and  described  in  the  previous  section,  has  some 
important  limitations:  it  is  applicable  only  to  robotic  networks  with  synchronized  and  reliable  communication 
along  all  edges  of  the  Delaunay  graph  (computed  as  a  function  of  the  robots’  positions).  In  other  words,  the 
law  ©  requires  that  there  exists  a  predetermined  common  communication  schedule  for  all  robots  and,  at  each 
communication  round,  each  robot  must  simultaneously  and  reliably  communicate  its  position.  Note  that  the 
Delaunay  graph,  interpreted  as  a  communication  graph,  has  the  following  drawbacks:  for  worst-case  robots’ 
positions,  a  robot  might  have  N  —  1  neighbors  in  the  Delaunay  graph  and/or  might  have  a  neighbor  that  is 
arbitrarily  far  inside  the  environment.  Therefore,  each  robot  must  be  capable  to  communicate  potentially  to  all 
other  robots  and/or  to  robots  at  large  distances. 

Given  this  broad  range  of  undesirable  limitations,  the  aim  of  this  paper  is  to  reduce  the  communication 
requirements  of  distributed  coverage  algorithms,  in  terms  of  reliability,  synchronization  and  topology.  Here  are 
some  relevant  questions  that  constitute  our  informal  problem  statement: 

Is  it  possible  to  optimize  robots  positions  and  environment  partition  with  asynchronous,  unreliable, 
and  delayed  communication?  Specifically,  what  if  the  communication  model  is  that  of  gossiping 
agents ,  that  is,  a  model  in  which  only  a  pair  of  robots  can  communicate  at  any  time?  Since  Voronoi 
partitions  generated  by  gossiping  and  moving  agents  cannot  be  computed  by  gossiping  agents,  how 
do  we  update  the  environment  partition? 

To  answer  these  questions,  the  next  subsections  propose  an  innovative  partition-based  gossip  approach,  in  which 
the  robots’  positions  essentially  play  no  role  and  where  instead  dominance  regions  are  iteratively  updated. 
Designing  coverage  algorithms  as  dynamical  systems  on  the  space  of  partitions  has  the  key  advantage  that  one 
is  not  restricted  to  working  only  with  Voronoi  or  anyway  position-dependent  partitions. 

Example  3.1  (The  Lloyd  algorithm  in  the  partition-based  approach)  The  distributed  coverage  algorithm  © 
updates  the  robots’  positions  so  as  to  incrementally  minimize  the  function  Tt\ oronoi,  while  the  environment  par¬ 
tition  is  a  function  of  the  robots  ’  positions.  In  this  paper  we  take  a  dual  approach:  we  consider  an  algorithm 
that  evolves  partitions.  From  this  partition-based  viewpoint ,  the  coverage  algorithm  is  an  iterated  map  on  Vn 
and  equation  ©  is  rewritten  as  v(t  +  1)  =  V(Cd(u(f))).  □ 

3.2  The  gossip  coverage  algorithm 

In  this  subsection  we  present  a  novel  partition-based  coverage  algorithm  in  which,  at  each  communication  round, 
only  two  regions  communicate.  Recall  the  notion  of  bisector  half-space  from  equation  ©. 
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Gossip  Coverage  Algorithm 

For  all  t  £  Z>o,  each  agent  i  £  {1,...,IV}  maintains  in  memory  a  dominance  region  Vi(t).  The  collection 
('Ui(O), . . . ,  wjv(O))  is  an  arbitrary  polygonal  IV-partition  of  Q.  At  each  t  £  Z>o  a  pair  of  communicating  regions, 
say  Vi(t)  and  Vj(t),  is  selected  by  a  deterministic  or  stochastic  process  to  be  determined.  Every  agent  k  £  {i,j} 
sets  Vk{t  +  1)  :=  Vk{t).  Agents  i  and  j  perform  the  following  tasks: 

1:  agent  i  transmits  to  agent  j  its  dominance  region  Vi(t)  and  vice-versa 
2:  both  agents  compute  the  centroids  Cd (vi(t))  and  Cd {vj(t)) 

3:  if  Cd(ui(t))  =  Cd (vj(t))  then 

4:  Vi(t  +  1)  :=  Vi(t)  and  Vj(t  +  1)  :=  Vj(t) 

5:  else 

6:  Vi(t  +  1)  :=  (vi(t)UVj(t))  n  ffbisectoI(Cd(vi(t));  Cd(vj (t))) 

Vj(t+  1)  :=  (vi(t)UVj(t))  nHUsectoIlGd(vj(t));Cd(vi(t))) 

In  other  words,  when  two  agents  with  distinct  centroids  communicate,  their  dominance  regions  evolve  as 
follows:  the  union  of  the  two  dominance  regions  is  divided  into  two  new  dominance  regions  by  the  hyperplane 
bisecting  the  segment  between  the  two  centroids;  see  Fig.  [3]  As  a  consequence,  if  the  centroids  Cd(uj(t)), 
Cd are  distinct,  then  {vi(t  +  1),  Vj(t+  1)}  is  the  Voronoi  partition  of  the  set  Vi(t)  U  Vj(t)  generated  by  the 
centroids  Cd(«j(i))  and  Cd(«j(t)). 


Figure  3:  The  gossip  coverage  algorithm.  The  left  and  right  figure  contain  the  initial  partition  and  the  partition 
after  one  application  of  the  gossip  coverage  algorithm.  In  the  middle  figure  we  show  the  two  centroids  and  (with 
a  dashed  line)  the  segment  determining  the  bisector  half-space. 


We  claim  that  the  algorithm  is  well-posed  in  the  sense  that  the  sequence  of  collections  {v{t)}tex>0  generated 
by  the  algorithm  is  an  IV-partition  at  all  times  t.  Indeed,  it  is  immediate  to  see  that  the  first  two  properties  in 
Definition  12. II are  satisfied  at  all  time  if  they  are  satisfied  at  initial  time.  Finally,  at  all  times  t,  each  component 
of  v(t)  is  clearly  closed  and  has  non-empty  interior  because  of  the  following  geometric  fact:  there  exists  no 
half-plane  containing  the  interior  of  a  region  and  not  containing  the  centroid  of  the  same  region. 

Now,  for  any  i,j  £  {1, . . . ,  N}  with  i  ^  j,  define  the  map  Tij  :  Vn  — ►  Vn  by 


Tij(v) 


v,  if  Cd(vi)  =  Cd(uj), 

(iq, . . .  ,Vi, . . .  ,Vj, . . .  ,vN),  otherwise, 


where 

Vi  —  (  c,;  U  Vj  'j  n  ^bisector  (Cd(^);  Cd(vj))  , 

Vj  =  (vi  LJ  Vj  j  n  //bisector  (Cd ('Cj ) ;  Cd('C7;)). 

The  dynamical  system  on  the  space  of  partitions  is  therefore  described  by,  for  t  £  Z>o, 


u(t-t-l)  =Tij(v(t)), 


(10) 


(11) 


8 


together  with  a  rule  describing  what  pair  of  regions  (i,j)  is  selected  at  each  time.  We  also  define  the  set-valued 
map  T  :  Vn  =4  Vn  by  T(v)  =  {Tiffv)  \  i,  j  £  {1, . . . ,  N},  i  ^  j}.  The  map  T  describes  one  iteration  of  the  gossip 
coverage  algorithm;  an  evolution  of  the  gossip  coverage  algorithm  is  one  of  the  solutions  to  the  non-deterministic 
set- valued  dynamical  system  v(t  +  1)  £  T(v(t)). 

Remark  3.2  (Gossip  disk-covering  control)  We  believe  that  our  gossip  and  partition-based  algorithmic  ap¬ 
proach  is  applicable  to  a  broad  range  of  coverage  problems.  For  example ,  the  worst-case  multicenter  function  m 
is  defined  by  H WOrst  {p,  v)  =  maxi6{i maxgei)i  ||<?  —  p,  ||.  Maximizing  HWOISt  is  equivalent  to  covering  Q  with 
N  disks  of  smallest  radius  centered  atp.  As  in  Proposition  for  any  v  £  Vn  and  p  £  Qn\Sn,  one  can  prove 
Tfworst (V(p) , p)  <  'Hv,orst{v,p)  and  HWorst(v, CC(u))  <  HWorSt(^,p),  where  CC(v)  is  the  array  of  circumcenters 
of  the  components  of  v.  Hence,  a  gossip  coverage  algorithm  for  7iworst  is  designed  by  replacing  centroid  with 
circumcenter  operations  in  m-  We  leave  this  and  further  extensions  to  future  works.  □ 


3.3  Analysis  results  for  the  gossip  coverage  algorithm 

In  this  subsection  we  state  the  main  analysis  and  convergence  results  for  the  gossip  coverage  algorithm. 

We  begin  by  studying  the  fixed  points  of  T  and  by  introducing  an  appropriate  cost  function  with  monotonicity 
properties  along  T.  Regarding  the  algorithm’s  fixed  points,  we  extend  Definition  12.51  as  follows.  A  partition  v 
is  mixed  centroidal  Voronoi  if,  for  all  pairs  ( Vi,Vj )  with  i  ^  j ,  either  Cd(ui)  =  Cd(uj)  or  ( Vi,Vj )  is  a  centroidal 
Voronoi  partition  of  Vi  U  Vj,  that  is,  Vi  =  ( Vi  Uu,)  H Rbisector(Cd(uj);  Cd(uj)). 

Lemma  3.3  (Fixed  points  of  T  and  centroidal  Voronoi  partitions)  For  i,j  £  {1, . . . ,  N},  j  ^  i,  denote 
the  set  of  fixed  points  ofTij  :  Vn  — >  Vn  by  Fij  =  {u  £  Vn  |  Tij(v )  =  v}.  The  following  statements  hold: 

(i)  rij^j  F^  equals  the  set  of  mixed  centroidal  Voronoi  partitions;  and 

(ii)  if  v  is  a  mixed  centroidal  Voronoi  partition  satisfying  Cd(uj)  ^  Cd(uj)  for  j  ^  i,  then  v  is  centroidal 
Voronoi. 


Next,  we  define  the  partition-based  multicenter  function  7Ycentroid  :  Vn  — >  R>o  by 

N 


^fcentroid(^)  —  multicenter  ('C ;  ^d(u))  —  ^  ^  j  /(||q  Cd(Uj  )  ||  ^)(j){(f)dq. 

i=  1  Jv* 


(12) 


Lemma  3.4  (Monotonicity  of  7iCentroid  along  T )  For  i,j  £  (1, . . . ,  N},  i  ^  j, 

centroid (u))  A  7YCentroid (’c) ;  far  all  V  £  VjV>  and 

^centroid {Tij{v))  <  Wcentroid (v) ,  iff  Tij{v)  and  v  differ  by  a  set  of  measure  zero. 

The  proofs  of  Lemmas  13 . 31  and  13 .41  consist  of  elementary  manipulations  and  are  omitted  in  the  interest  of  brevity. 
In  short,  we  have  established  that  the  function  7YCentroid  monotonically  decreases  along  each  Ty  when  away  from 
fixed  points,  and  that  centroidal  Voronoi  partitions  are  the  fixed  points  of  all  Tij  provided  centroids  are  distinct. 

We  now  prepare  to  state  the  main  convergence  result  for  T.  We  need  to  introduce  some  useful  properties 
for  sequences  of  partitions  and  for  switching  signals. 

Definition  3.5  (Non-degeneracy)  A  sequence  of  N -partitions  {v(t)}t^z>0  is 

(i)  (uniformly)  distinct  centroidal  if  there  exists  e  >  0  such  that,  for  all  t  £  Z>o  and  iff  £  {1, . . . ,  N},  i  j , 

one  has  ||  Cd(uj(t))  —  Cd(uj(t))||  >  e; 

(ii)  (uniformly  componentwise)  non- vanishing  if  there  exists  e  >  0  such  that,  for  all  t  £  Z>o  and  i  £  {1, . . . ,  N}, 
the  Lebesgue  measure  of  Vi(t)  is  greater  than  e;  and 

(Hi)  (uniformly  componentwise)  finitely  convex  if  there  exists  t  £  N  such  that,  for  all  t  £  Z>o  and  i  £ 
{1, . . . ,  N},  the  set  Vi(t)  is  the  union  of  at  most  £  convex  sets. 
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Moreover,  the  sequence  v  is  said  to  be  non-degenerate  if  it  is  distinct  centroidal,  non-vanishing  and  finitely 
convex. 

For  example,  a  sequence  of  partitions  is  finitely-convex  if  each  component  of  each  partition  in  the  sequence 
is  the  union  of  a  uniformly  bounded  number  of  polygons  with  a  uniformly  bounded  number  of  vertices. 

Definition  3.6  (Uniform  and  random  persistency)  Let  X  be  a  finite  set. 

(i)  A  map  a  :  Z>o  — >  X  is  uniformly  persistent  if  there  exists  a  duration  A  £  N  such  that,  for  each  x  £  X , 
there  exists  an  increasing  sequence  of  times  {tk}kez>0  C  Z>o  satisfying  tk+i  —  tk  <  A  and  a(tk)  =  x  for 
all  k  £  Z>o- 

(ii)  A  stochastic  process  a  :  Z>o  —>  X  is  randomly  persistent  if  there  exists  a  probability  p  £  ]0,1[  and  a 
duration  A  £  N  such  that,  for  each  x  £  X  and  for  each  t  £  Z>o,  there  exists  k  £  (1, . . . ,  A}  satisfying 

P[a(t  +  k)  =  x  |  aft), . . . ,  cr(l)]  >  p. 

We  are  now  ready  to  state  the  main  deterministic  and  stochastic  convergence  results  for  the  gossip  coverage 
algorithm.  It  is  convenient  to  postpone  to  Section  16.21  the  theorem  proof  and  the  definition  of  convergence  in 
the  space  of  partitions. 

Theorem  3.7  (Convergence  under  persistent  gossip)  Consider  the  gossip  coverage  algorithm  T  and  the 
evolutions  {v(t)}t^i>0  C  Vn  defined  by 

v(t+l)  =  Ta{t)(v(t)),  for  t  £  Z>o, 

where  a  :  Z>o  — >  {(i,j)  £  {1, . . . ,  TV}2  |  i  ^  j}  is  either  a  deterministic  map  or  a  stochastic  process.  Then  the 
following  statements  hold: 

(i)  if  a  is  a  uniformly  persistent  map,  then  each  non-degenerate  evolution  v  converges  to  the  set  of  centroidal 
Voronoi  partitions;  and 

(ii)  if  a  is  a  randomly  persistent  stochastic  process,  then  each  evolution  v,  conditioned  upon  being  non¬ 
degenerate,  converges  almost  surely  to  the  set  of  centroidal  Voronoi  partitions. 

The  statements  of  Theorem  13. 71  relv  upon  the  assumption  of  non-degenerate  evolutions.  It  is  our  conjecture 
that,  starting  from  generic  polygonal  partition,  this  assumption  is  typically  satisfied.  We  will  discuss  some 
numerical  evidence  to  this  effect  in  the  next  section. 

Lemma  13.41  indicates  how  the  function  TUentroid  plays  the  role  of  a  Lyapunov  function  for  the  dynamical 
system  defined  by  T.  To  provide  a  complete  Lyapunov  convergence  proof  of  Theoreml3.7l  we  set  out  to  establish 
three  sets  of  relevant  results.  First,  we  need  to  establish  extensions  of  the  Krasovskii-LaSalle  invariance  principle 
for  set-valued  dynamical  systems  over  compact  metric  spaces.  Second,  we  need  to  establish  the  compactness 
properties  of  the  space  of  non-degenerate  partitions.  Third,  we  need  to  establish  the  continuity  of  the  relevant 
geometric  maps.  These  three  topics  are  the  subjects  of  Section  [4]  [5]  and  [0  respectively. 

3.4  Designing  an  adjacency-based  and  continuous  algorithm 

The  gossip  coverage  map  T  has  the  undesirable  feature  that  it  entails  communication  exchanges  between  any 
two  regions.  This  communication  requirement  might  be  too  onerous  in  some  multiagent  applications.  Ideally  we 
would  like  to  require  communications  only  between  adjacent  regions,  that  is,  regions  whose  boundaries  touch, 
or  between  “nearby”  regions.  We  believe  such  communication  requirements  may  be  easily  achieved  in  robotic 
networks  and  we  detail  a  sample  implementation  for  robots  with  range-dependent  Poisson  communication  in 
Appendix  [A]  Additionally,  we  aim  to  design  a  continuous  gossip  coverage  map.  We  require  the  modified  map 
to  be  continuous  for  technical  reasons:  the  invariance  principles  we  adopt  for  the  convergence  analysis  require 
continuity  of  the  dynamical  system. 
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Motivated  by  this  discussion,  we  modify  the  map  T  to  rely  upon  only  adjacency-based  communication  and 
to  be  continuous.  First,  we  introduce  a  pseudodistance  notion  between  sets.  Given  closed  A  C  Q  and  B  C  Q 
with  non-empty  interior,  define 


pseudodist(A,  B)  =  inf{||a  —  6||  |  (a,  b)  £  int  (A)  x  int(.B)}. 

Second,  we  select  a  positive  constant  <5  <C  diam(Q)  and  denote  by  Ts  :  Vn  —■ •  Vat  the  modified  gossip  coverage 
map  to  be  defined  in  what  follows.  For  any  i,j  £  {1, . . . ,  TV},  i  y  j,  we  give  the  following  partial  definition: 


v,  if  (||  Cd(«<)  -  Cd^OII  =  o)  or  ( pseudodist(r>j,  Vj)  >  S), 

T{v),  if  (||  Cd(ui)  —  Cd(uj)||  >  <5)  and  ( pseudodist(i>i,  Vj)  =  0). 


(13) 


Therefore,  if  either  Cd(?7)  and  Cd(uj)  coincide  or  the  pseudodistance  between  Vi  and  Vj  is  larger  than  <5,  then 
Tf-  (v)  =  v,  that  is,  the  map  leaves  the  partition  unchanged.  Additionally,  if  the  pseudodistance  between  the 
regions  vl  and  Vj  is  zero  (vi  and  Vj  are  adjacent)  and  the  distance  between  Cd(i>i)  and  Cd(t!j)  is  larger  than  <5, 
then  Tfj(v)  =  Tij(v). 

Next,  we  consider  partitions  that  do  not  satisfy  either  of  the  two  conditions  in  definition  (1131).  We  define  the 
unit  saturation  function  sat  :  R>o  — i ►  [0, 1]  by  sat(x)  =  x  if  x  £  [0,1],  and  sat(x)  =  1  if  x  >  1  and  the  scaling 
function  0ij  :  Vn  — ►  [0, 1]  by 

i %j{v )  =  sat ( 1 1  Cd(uj)  —  Cd(i>j)||/<5)  ^1  —  sat(pseudodist(vj,  Vj)/S)^ . 

The  first  condition  and  the  second  condition  in  (1131)  correspond  precisely  to  (3 ij(v)  =  0  and  (3ij{v)  =  1,  respec¬ 
tively.  For  partitions  v  satisfying  0  <  Pij(v)  <  1,  we  aim  to  define  Ts  so  as  to  continuously  interpolate  between 
the  identity  map  and  the  map  T;  see  Fig.  0]  for  an  illustration.  Let  Ri  =  Vi  D  l?bisector(Cd(uj),  Cd(uj))  and 


Figure  4:  Modified  gossip  between  close  but  not  adjacent  regions  (0  <  (3 12 (u)  <  1).  The  bisecting  line  7_l 
borders  the  set  R2  =  V2  D  -ffbisector(Cd(?;i),  Cd(i>2))  that,  in  the  map  T,  is  assigned  to  v\  (see  Fig.  [3]).  According 
to  Ts  instead,  only  the  subset  R2  C  R2  is  assigned  to  v\.  Loosely  speaking,  the  “width”  of  R2  equals  /3i2(v) 
times  the  “width”  of  ffe,  where  “width”  of  a  set  is  the  maximum  distance  from  a  point  in  the  set  to  7_l. 


Rj  =  Vj  n7Lbisector(Cd(rii),  Cd(uj)).  Define  the  line  7_l  =  9i?bisector(Cd(uj),  Cd(ui))  and 

Pi  =  a  point  in  int  {Ri)  that  is  maximally  distant  from  7_l,  ^ 

Pj  =  a  point  in  int  (Rj)  that  is  maximally  distant  from  7j_. 

Next,  note  that  for  each  q  £  Ri  U  Rj  there  exists  a  unique  line,  say  7g,  that  is  parallel  to  7_l  and  passes  through 
q.  Based  on  this  notion,  we  define 

Ri  =  {q  £  Ri  I  dist(pj,  7g)  <  (3ij{v)  dist(pi,  7x)  or  dist(g,7_i_)  >  dist(pj,  7_l)}, 

Rj  =  {q  £  Rj  I  dist (j>j,jg)  <  /3lj(u)dist(pi,7_L)  or  dist(g,7_i_)  >  dist(pj,7j_)}. 

We  can  now  complete  the  partial  definition  (|  1 3 1) .  For  all  v  with  0  <  /3ij(v)  <  1,  that  is,  for  all  partitions  not 
already  dealt  with  in  definition  (I13|) ,  we  define 

T-j{v)  =  (vi, . . . ,  (vi  \  Ri)  U  Rj, . . . ,  (vj  \  Rj)  U  vN). 

^  -  -  s  ^  s 

■  V  ■  ■  V  "  ^ 

ith  entry  j th  entry 
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As  discussed  for  T,  one  can  prove  that  the  map  Ts  :  Vn  =4  Vn  defined  by  Ts{v )  =  {T^(v)  \  i,  j  G  {1, . . . ,  N},  i  ^ 
j},  is  well-posed  and  has  the  following  properties. 

Theorem  3.8  (Convergence  of  modified  gossip  map)  Consider  the  modified  gossip  coverage  algorithm  Ts 
and  the  evolutions  {v(t)}t^z>0  C  Vn  defined  by 

v(t+  1)  =  for  t  G  Z>o, 

where  a  :  Z>o  — ►  {(*,_))  G  { 1, . . . ,  iV}2  |  i  ^  j}  is  either  a  deterministic  map  or  a  stochastic  process.  Then  the 
following  statements  hold: 

(i)  if  cr  is  a  uniformly  persistent  map ,  then  each  non-vanishing  and  finitely- convex  evolution  v  converges  to 
the  set  of  mixed  centroidal  Voronoi  partitions;  and 

(ii)  if  a  is  a  randomly  persistent  stochastic  process,  then  each  evolution  v,  conditioned  upon  being  non-vanishing 
and  finitely  convex,  converges  almost  surely  to  the  set  of  mixed  centroidal  Voronoi  partitions. 

3.5  Simulation  results  and  implementation  remarks 

We  have  extensively  simulated  the  partition-based  gossip  coverage  algorithm  T  on  a  2-dimensional  polygonal 
environment  with  uniform  density  and  performance  function  f{x)  =  x2 .  Simulations  have  been  implemented 
as  a  Matlab  program,  using  the  General  Polygon  Clipper  Library  to  perform  operations  on  polygons.  We 
adopted  the  following  communication  model:  at  each  iteration,  a  region  pair  is  chosen,  uniformly  at  random, 
among  all  pairs  of  adjacent  regions.  Fig.  [5]  is  an  illustration  of  a  typical  evolution. 


Figure  5:  An  example  simulation  of  the  gossip  coverage  algorithm  with  uniform  random  edge  selection.  The 
environment  Q  is  a  rectangle  with  uniform  density,  centroids  are  computed  with  performance  function  f(x)  =  x2, 
and  N  =  6  regions  compose  the  partition.  Snapshots  of  an  evolution  are  shown  for  t  G  {0, 20,  50, 100,  300}.  One 
may  verify  numerically  that  the  sequence  converges  asymptotically  to  a  centroidal  Voronoi  partition.  At  t  =  20 
one  of  the  regions  is  disconnected. 

Our  first  numerical  finding  is  that  the  gossip  coverage  algorithm  appears  to  converge  to  centroidal  Voronoi 
partitions  from  all  initial  conditions.  This  is  the  same  property  that  the  Lloyd  synchronous  coverage  algo¬ 
rithm  ©  is  known  to  possess.  In  other  words,  our  numerically-computed  sequences  of  partitions  always  converge 
to  centroidal  Voronoi  partitions  -  even  though  our  theoretical  analysis  (1)  requires  a  continuous  interpolation 
from  T  to  Ts  and  (2)  does  not  exclude  convergence  to  degenerate  partitions  where  some  component  regions 
might  have  coincident  centroids,  or  might  have  empty  interiors,  or  might  be  composed  of  “polygons  with  an 
infinite  number  of  vertices,”  that  is,  arbitrary  sets. 

A  second  numerical  finding  is  that,  throughout  numerous  sample  executions,  the  resulting  polygonal  regions 
rarely  have  complicated  shapes  and  large  numbers  of  vertices.  This  is  good  news  because  of  our  assumption 
of  finite  convexity  and  because  large  numbers  of  vertices  affect  both  the  computation  and  the  communication 
burden  of  the  gossip  coverage  algorithm. 

Finally,  it  is  possible,  and  we  have  observed  it  numerically,  to  have  evolutions  of  the  algorithm  that,  before 
converging  to  centroidal  Voronoi  partitions,  have  components  with  disconnected  regions.  We  believe  that  there 
might  be  applications  where  it  is  desirable  to  maintain  connectivity  of  the  components  of  the  partition  and, 
therefore,  we  sketch  here  how  to  design  a  connectivity-preserving  algorithm.  Note  that  the  update  step  of  the 
partition-based  coverage  algorithm  amounts  to  the  exchange,  among  the  agents,  of  a  region,  which  consists 
in  general  of  several  connected  components.  In  the  connectivity-preserving  algorithm,  such  components  are 
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considered  individually,  and  each  of  them  is  traded  only  if  this  can  be  done  without  loosing  connectivity;  if  not, 
the  component  is  kept  by  the  previous  owner.  Numerical  simulations  indicate  that  such  an  algorithm  leads  to 
centroidal  Voronoi  partitions  as  well. 

4  On  the  Krasovskii— LaSalle  invariance  principle:  set-valued  maps 
on  metric  spaces 

In  this  section  we  consider  discrete-time  set-valued  dynamical  system  defined  on  metric  spaces.  Our  goal  is  to 
provide  some  extensions  of  the  classical  Krasovskii-LaSalle  Invariance  Principle;  we  refer  the  interested  reader 
to  0  ffiJl  134]  for  recent  invariance  principles  for  switched  continuous-time  and  hybrid  systems. 

We  start  by  reviewing  some  preliminary  notions  including  set-valued  dynamical  systems,  continuity  and 
invariance  properties,  and  Lyapunov  functions.  On  a  metric  space  ( X ,  d ),  where  X  is  a  set  and  d  is  a  metric  on 
X ,  a  set- valued  map  T  :  X  X  is  non-empty  if  T(x )  ^  0  for  all  x  £  X.  An  evolution  of  the  dynamical  system 
determined  by  a  non-empty  set-valued  map  T  is  a  sequence  {xn}n£z>0  C  X  with  the  property  that 

Xn-\-l  £  T(xn ),  Tl  £  Z>0- 

In  other  words,  we  regard  a  set-valued  map  as  a  nondeterministic  discrete-time  dynamical  system.  For  set-valued 
maps  we  introduce  notions  of  continuity  and  invariance  as  follows.  A  set-valued  map  T  is  closed  at  x  £  X  if, 
for  all  pairs  of  convergent  sequences  Xk  — >  x  and  x'k  — ►  x'  such  that  x'k  £  T{xk ),  one  has  that  x'  £  T(x). 

Additionally,  T  is  closed  on  W  C  X  if  it  is  closed  at  all  w  £  W.  A  set  W  C  X  is  weakly  positively  invariant  for 
T  if  T{w )  n  W  is  non-empty  for  all  w  £  W.  A  set  W  is  strongly  positively  invariant  for  T  if  T(w)  C  W  for  all 
w  £  W. 

We  are  ready  now  to  state  a  Krasovskii-LaSalle  invariance  principle  for  set-valued  maps  defined  on  metric 
spaces.  This  result  extends  the  Global  Convergence  Theorem  in  [35]  to  more  general  Lyapunov  functions.  Its 
proof  follows  the  lines  of  the  proof  of  Theorem  1.21  in  [5],  and  is  thus  omitted. 

Lemma  4.1  (Krasovskii-LaSalle  invariance  principle  for  set-valued  maps)  Let(X,d)  be  a  metric  space 
and  T  :  X  z4  X  be  non-empty.  Assume  that: 

(i)  there  exists  a  compact  set  W  C  X  that  is  strongly  positively  invariant  for  T; 

(ii)  there  exists  a  function  U  :  W  — >  R  such  that  U{w')  <  U(w),  for  all  w  £  W  and  w'  £  T(w); 

(in)  the  function  U  is  continuous  on  W  and  the  map  T  is  closed  on  W. 

Then  there  exists  c  £  R  such  that  each  evolution  of  T  with  initial  condition  in  W  approaches  a  set  of  the  form 
S'nt/_1(c),  where  S  is  the  largest  weakly  positively  invariant  set  contained  in 

{w  £W\3w'  £  T(w)  such  that  U(w')  =  U(w)}. 

In  this  paper,  given  the  metric  space  ( X ,  d),  we  deal  with  a  set- valued  map  T  :  X  X  defined  by  a  collection 
of  maps  Ti, ... ,  Tm  :  X  — >■  X  via  the  equality  T(x)  =  {Ti(x), . . . ,  Tm{x )}  for  x  £  X .  For  this  kind  of  set-valued 
maps,  closedness  is  related  to  continuity  of  ordinary  maps. 

Lemma  4.2  Let  T\1 . . .  ,Tm  :  X  — >  X  be  continuous  on  W  C  X.  The  set-valued  map  T  :  X  X  defined  by 
T(x)  =  {Ti(x), . . . ,  Tm(x)}  is  closed  on  W. 

Proof.  Let  wn  — >  w  and  w'n  — >  w'  be  a  pair  of  convergent  sequences  in  W,  such  that  w'n  £  T(wn).  We  claim 
that  the  continuity  of  Ti, . . . ,  Tm  implies  w'  £  T(w). 

Note  that,  by  hypothesis,  for  all  n  £  Z>o  there  exists  in  £  (1, . . .  ,m)  such  that  w'n  =  Tin(wn).  Because  the 
set  {1, . . . ,  to}  is  finite,  there  exists  an  index  j  £  {1, . . . ,  to}  that  appears  infinitely  many  times  in  the  sequence 
{*n}n£N-  Consider  the  subsequences  {wni}  C  {w„}  and  {ui'ni}  C  {w'n},  such  that  w'ni  =  Tj(wni).  Clearly,  we 
have  that  wni  — ►  w  and  w'  — >  w' ,  where  from  the  continuity  of  T)  it  follows  that  w'  =  Tj(w).  Thus,  w'  £  T(w) 
and  the  claim  is  proved.  I 

The  following  result  is  a  stronger  version  of  Lemmal4~Tl  for  a  particular  class  of  set- valued  dynamical  systems. 
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Theorem  4.3  (Uniformly  persistent  switches  imply  convergence)  Let  ( X ,  d)  be  a  metric  space.  Given 
a  collection  of  maps  Ti, . . . ,  Tm  :  X  — ►  X,  define  the  set-valued  map  T  :  X  =1  X  by  T( x)  =  {T\(x), . . . ,  Tm(x)} 
and  let  {xn}n^z>0  be  an  evolution  ofT.  Assume  that: 

(i)  there  exists  a  compact  set  W  C  X  that  is  strongly  positively  invariant  for  T; 

(ii)  there  exists  a  function  U  :  W  — >  R  such  that  U(w')  <  U(w),  for  all  w  £  W  and  w’  £  T(w)  \  {w}; 

(in)  the  maps  Ti,  for  i  £  {1, . . . ,  to},  and  U  are  continuous  on  W;  and 

(iv)  for  all  i  £  there  exists  an  increasing  sequence  of  times  {nk\  k  £  Z>o}  such  that  xnk+i  = 

Ti(xnk )  and  (nk+i  —  Uk)  is  bounded. 

If  Xq  £  W,  there  exists  cel  such  that  the  evolution  {xn}nez>0  approaches  the  set 

(F1n---nFm)nt/-1(c), 

where  Fi  =  {w  £  W  \  Ti(ui)  =  u>}  is  the  set  of  fixed  points  of  Ti  in  W,  i  £  (1, . . . ,  to}. 

Loosely  speaking,  (i)  the  compactness  of  a  strongly  positively  invariant  set,  (ii)  a  monotonicity  property  for  a 
Lyapunov  function,  (iii)  continuity  properties,  and  (iv)  uniformly  persistent  switches  among  finitely  many  maps, 
together  ensure  convergence  of  each  evolution  to  the  intersection  of  the  fixed  points  of  the  maps. 

Proof.  [Proof  of  Theorem  14. 3j  Let  S  be  the  largest  weakly  positively  invariant  set  contained  in 

{w  £  W  |  3  w'  £  T{w)  such  that  U(w')  =  U(w)}  =  Fi  U  •  •  ■  U  Fm. 

Since  T  is  closed  by  Lemma  14.21  the  assumptions  of  Lemma  14.11  are  met;  hence  there  exists  c  £  R  such  that 
U(xn )  — ►  c  and  xn  — >  S  n  C/_1(c). 

Let  fl(xn)  denote  the  w-limit  set  of  the  sequence  { xn  \  n  £  Z  >0}.  If  we  show  that  f 2(x„)  C  (Fi  fl  •  •  •  PI  Fm)  fl  C/_1(c), 
then  the  statement  of  the  theorem  is  proved.  We  proceed  by  contradiction.  To  this  aim,  let  x  £  S'n[/_1(c)  \ 

((Fi  n  •  •  •  n  Fm)  n  t/_1(c))  and  let  {xnh}h£Z>0  be  a  subsequence  such  that  xnh  —>  x. 

Observe  that  for  each  x  £  S  \  (Fi  fl  ■  •  •  f~lFm),  there  exists  a  non-empty  set  C  {1, . . . ,  to}  such  that, 
x  =  Ti  (ai)  if  i  £  Zg,  and  x  ^  Ti  ( x )  if  i  X By  the  continuity  of  the  maps  Ti ,  there  exists  <5  £  R>o  such  that, 
if  i  fi  X%,  then  Ti(x)  ^  x  for  all  x  £  Bs  (x)  =  {x  £  W  \  d(x,  x)  <  S}.  Let  now 

7 s  =  min  i  min  (u  (x)  —  U  (TAx))  )  1  >  0. 

(leBi(J)  V  /  J 

By  hypothesis,  if  %  f  T%,  then  U(Ti(x))  <  U(x)  for  all  x  £  B$(x).  Hence,  since  B$( x)  is  closed,  and  U  and  the 
maps  Ti  are  continuous,  we  deduce  that  7  >  0. 

Observe  now  that  hypothesis  (iv)  implies  the  existence  of  a  duration  D  £  N  such  that  every  map  Ti, 
i  £  {l,...,m},  is  applied  at  least  once  within  every  interval  [n,n  +  D[,  for  all  n  £  Z>o-  Consider  the  set 
{Ti}igi  ;  is  a  collection  of  continuous  maps  having  x  as  fixed  point.  Then,  there  exists  a  suitable  e  £  R>o 
such  that,  given  any  r-uple  (j  1, . . . ,  jr)  £  XI,  r  <  D,  we  have  that  Tj1  o  Tj2  o  Tj3  o  . . .  o  Tjr(x)  £  Bs(x)  for  all 
w  £  Be(x). 

Select  now  k  such  that  the  element  xnk  in  the  subsequence  {xnh  }hez>0  satisfies  d  ( xnk ,  x)  <  e  and  U ( xnk )  —  c  <  75. 
Let 

s  =  min{t  £  [1,  D[  \  3j  f  X3  such  that  xnk+t+i  =  Tj  (; xnk+t )}. 

Observe  that  U(xnk+S)  -  c  <  js  and  U(xnk+S )  -  U(Tj(xnk+s))  >  7 s  implying  that  U  (Tj  (xUk+s))  <  c.  This  is 
a  contradiction.  I 

Remark  4.4  An  alternate  proof  of  this  theorem  can  be  given  by  applying  an  invariance  principle  obtained  in  | 
on  an  appropriately- designed  dynamical  extension  of  the  discrete-time  set-valued  system.  □ 

In  Appendix[B]we  show  how  to  persistent  switching  assumption  is  necessary.  Next,  we  provide  a  probabilistic 
version  of  the  previous  theorem. 
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Theorem  4.5  (Persistent  random  switches  imply  convergence)  Let  (X,d)  be  a  metric  space.  Given  a 
collection  of  maps  Ti, . . . ,  Tm  :  X  — ►  X,  define  the  set-valued  map  T  :  X  =4  X  by  T(x)  =  {T\{x), . . . ,  Tm(x)}. 
Given  a  stochastic  process  a  :  Z>o  — »  {1, . . .  ,m},  consider  an  evolution  {xn}n&i>0  ofT  satisfying 

-^cr(n)  (^n  )  ■ 


Assume  that: 


(i)  there  exists  a  compact  set  W  C  X  that  is  strongly  positively  invariant  for  T ; 

(ii)  there  exists  a  function  U  :  W  — >  R  such  that  U(w')  <  U(w),  for  all  w  £  W  and  w'  £  T(w)  \  {1/;}; 

(in)  the  maps  Ti,  for  i  £  {1, . . . ,  to},  and  U  are  continuous  on  W ;  and 

(iv)  there  exists  p€]0,l[  and  k  £  N  such  that,  for  all  i  £  {1, ,  m}  and  n  £  Z>o,  there  exists  h  £  (1, . . . ,  k} 
such  that 

P [er(n  +  h)  =  i  \  cr(n), . . . ,  <r(l)]  >  p. 

If  xq  £  W,  then  there  exists  c?l  such  that  almost  surely  the  evolution  {xn}nez>0  approaches  the  set 

(Pin---nFm)nt/-1(c), 


where  Fi  =  {rt>  £  W  \  Ti(w)  =  w}  is  the  set  of  fixed  points  of  Ti  in  W,  i  £  (1, . . . ,  to}. 

Loosely  speaking,  (i)  the  compactness  of  a  strongly  positively  invariant  set,  (ii)  a  monotonicity  property  for  a 
Lyapunov  function,  (iii)  continuity  properties,  and  (iv)  persistent  random  switches  among  finitely  many  maps, 
together  ensure  convergence  of  each  evolution  to  the  intersection  of  the  fixed  points  of  the  maps. 

Proof.  [Proof  of  Theorem  14. 5j  If  xq  £  W,  then  the  stochastic  process  a  induces  a  stochastic  process  taking 
values  in  W.  From  now  on,  we  restrict  our  attention  to  sequences  {xn}nez>0  such  that  xq  £  W .  In  other  words 
we  assume  that  the  sample  space  containing  all  the  evolutions  of  our  interest  is  given  by 

A  =  {{xn}n£Z>0  I  Xn  £  W  for  all  n  £  Z>o}. 

Let  S  be  the  largest  weakly  positively  invariant  set  contained  in 


{w  £  W  |  3w'  £  T(w)  such  that  U(w')  =  U(w)}  =  Fi  U  •  •  •  U Fm. 

From  Lemma  ITT!  we  know  that  there  exists  c  £  R  such  that  xn  — >  S  H  C/_1(c).  This  implies  that  the  following 
event  is  certain: 

E  =  i{^ra}nez>n  |  3  c  £  K  such  that  lim  U(xn)  =  c}. 

Let  fl(xn)  denote  the  w-limit  set  of  the  sequence  {xn  \  n  £  Z  >0}.  If  we  show  that  fl(xn)  C  ((Fi  fl  •  •  •  D  Fm)  ft  U~l  (c)) 
almost  surely,  then  the  statement  of  the  theorem  is  proved.  Next,  consider  the  event 


Ei  =  {{x„}  nez>0  |  3  x  £  S  \  (Fi  n...n  Fm)  such  that  x  £  LI  (x„)  } . 

Assume  by  contradiction  that  P  [E\]  >0.  Now  we  compute  P  [F|Fi].  Note  that,  for  each  x  £  S\ (F\  H  ...  H  Fm), 
there  exists  a  non-empty  set  Xj  C  {1, . . . ,  m}  such  that,  x  =  Ti  (x)  if  i  £  X^,  and  x  ^  Ti  (x)  if  i  ^  X^.  Similarly 
to  the  proof  of  Theorem  14.31  we  can  associate  to  each  x  a  positive  real  number  S  such  that  the  inequality 
x  ^  T  (x)  holds  true  for  all  x  £  Bs(x)  =  {x  £  W  \  d(x,  x)  <  5}  and  for  all  i  ^  X^.  Moreover,  we  can  define 


7  s  =  mm 


x£Bs(x] 


(U(x)  -  U(Ti(x )))  ]  , 


where  the  continuity  of  the  maps  X) ,  j  £  {1, . . .  ,to},  and  U,  and  the  closedness  of  the  set  Bs(x)  ensure  that 
15  >  0. 

Consider  the  set  {X,}ig:2>;  this  is  a  collection  of  continuous  maps  having  x  as  fixed  point.  Therefore,  there 
exists  a  suitable  e  £  E>o  such  that,  given  any  r-uple  (j\ , . . . ,  jr)  £  XI,  r  <  fc,  we  have  Tjx  o Tj2  o Tj3  o. .  .0 Tjr  (x)  £ 
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Bg(x)  for  all  x  £  Be(x).  Given  {xn}ne%>0,  if  there  exists  x  £  S  \  (Fi  fl . . .  fl  Fm)  such  that  x  £  fl  (xn),  then 
there  must  exist  { rih\  h  £  Z>o}  such  that  xnh  £  Be(x)  for  all  h  £  Z>o-  Moreover,  without  loss  of  generality  we 
can  assume  that  rih+i  —  Uh  >  k  for  all  h  £  Z>o-  Consider  now  the  event 

E3  =  {{in  £  {1,. .  .,m}  |  n  £  Z>0}  | 

3  h  such  that  iUh+s  £  I for  all  s  £  {1, . . . ,  k}  and  h  >h{. 


To  compute  P[£a],  we  define,  for  j  £  Z>o, 

E3 j  =  {{*„  £  {1, . . .  ,  m}  |  n  £  Z>0}  |  inh+s  £  h  for  all  s  £  {1, ... ,  k}  and  h>  j}. 

Observe  that  {E3j  \  j  £  Z>o}  is  a  countable  collection  of  disjoint  sets  such  that  E3  =  [jJLo  E3,j •  By  hypothesis 
we  have  that 

i 

P[^3j]  <  lim  TT(1  ~p)  =  0, 

l — *00  A  A 
s=j 

and  therefore  Pj-Ey  =  0.  This  implies  that,  almost  surely,  there  exists  a  subsequence  {rihs} sez.>0  1=  {nh}he,z>0 
with  the  property  that,  for  all  s  £  Z>o,  xnh  +i  =  Tj( xnh  )  for  some  i  £  and,  therefore,  also  with  the 
property  that  U(xnhs)  —  U(xnhk+ 1)  >  7 5.  Consequently,  almost  surely,  we  have  that  lim^oo  U(xnhs)  =  —00 
thus  violating  the  fact  that  E  is  a  certain  event.  This  implies  that  P  [Ei]  =  0  and  that,  almost  surely,  xn  — > 
(Fin...nFm)nc/_1(c).  '  ■ 

Remark  4.6  The  assumption,  in  Lemma\4-.1\  and  Theorems\4-3\  and\4-5\  that  W  is  strongly  positively  invariant 
ensures  that  any  evolution  of  T  with  initial  condition  in  W  remains  in  W.  By  relaxing  this  assumption,  it 
is  possible  to  provide  weaker  versions  of  these  results.  Specifically,  requiring  W  to  be  only  compact  and  not 
necessarily  strongly  positively  invariant,  the  thesis  of  Lemma  \4-l\  and  Theorems  \ 4-3\  and  |^.5|  do  not  hold  in 
general  for  any  evolution  of  T  with  initial  condition  on  W ,  but  are  still  valid  for  those  evolutions  {xn}n£z>0 
that  take  values  in  W  for  all  n  £  Z>o-  □ 

5  On  the  topology  of  the  space  of  partitions:  compactness  properties 
in  the  symmetric  difference  metric 

Motivated  by  the  invariance  principles  presented  in  Section[4l  we  study  metric  structures  on  the  set  of  partitions, 
with  a  focus  on  compactness  and  continuity  properties.  Specifically,  we  show  how  a  particular  subset  of  the 
set  of  partitions  can  be  regarded  as  a  compact  metric  space  and  how  certain  relevant  maps  are  continuous  over 
that  subspace.  In  this  section,  and  only  in  this  section,  the  assumptions  on  Q  are  relaxed  to  give  more  general 
results:  we  assume  that  Q  C  is  compact  and  connected  and  has  non-empty  interior. 

Let  C  denote  the  set  of  closed  subsets  of  Q.  We  would  like  to  introduce  a  topology  on  C  with  two 
properties:  C  is  compact  and  the  Voronoi  map,  the  centroid  map,  and  the  multicenter  function,  defined  in 
equations  ©  ©,  and  m  respectively,  are  continuous  over  C.  A  natural  candidate  is  the  topology  in¬ 
duced  by  the  well-known  jSJj  Hausdorff  metric  on  C :  given  two  sets  A,  B  £  C,  their  Hausdorff  distance  is 
dn{A,B)  =  max  {maxQg^  rnin^gs  d(a,  6),  max^gs  minagj4  d(a,  fe)}.  This  metric  induces  a  topology  on  C  which 
makes  it  a  compact  space,  but  is  not  suitable  for  our  purpose  because,  with  respect  to  this  topology,  the  Voronoi 
map,  the  centroid  map,  and  the  multicenter  function,  are  not  continuous;  see  Appendix  [Cl  Additionally,  note 
that,  unlike  the  Hausdorff  metric,  the  centroid  map  and  the  multicenter  function  are  insensitive  to  sets  of 
measure  zero. 

In  what  follows,  we  introduce  the  symmetric  difference  metric,  as  a  metric  insensitive  to  sets  of  measure  zero. 
Given  two  subsets  A,  B  £  C,  we  define  their  symmetric  difference  by  AqB  =  (AuB)\(A  fl  B).  Moreover,  letting 
p  denote  the  Lebesgue  measure  on  we  define  the  symmetric  difference  distance,  also  called  the  symmetric 
distance  for  simplicity,  c?q  :  C  x  C  ®>o  by 


de{A,  B)  —  p{A  0  B), 
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that  is,  the  symmetric  distance  between  two  sets  is  the  measure  of  their  symmetric  difference.  Given  these 
notions,  it  is  useful  to  identify  sets  that  differ  by  a  set  of  measure  zero:  we  write  A  ~  B  whenever  p(AQB)  =  0. 
Clearly,  ~  is  an  equivalence  relationship  on  C  and,  accordingly,  we  let  C*  =  C /  ~  denote  the  quotient  set  of 
closed  subsets  of  Q.  Now,  for  any  two  elements  A*  and  B*  of  C* ,  we  define  dg(A*,I?*)  =  dg(A,  B)  where  A 
and  B  are  any  representatives  of  A*  and  B* ,  respectively.  With  this  notion  of  dg  on  C*  x  C*,  it  is  easy  to  verify 
that  (C*,dg)  is  a  metric  space.  However,  to  the  best  of  our  knowledge  no  compactness  result  is  available  for 
(C*,de). 

Next,  we  introduce  a  particular  family  of  subsets  of  C  whose  structure  is  sufficiently  rich  for  our  algorithm. 
For  leN,  let  C  C  denote  the  set  of  £-conve x  and  closed  subsets  of  Q ,  that  is,  the  set  of  subsets  of  Q  equal 
to  the  union  of  l  convex  and  closed  subsets  of  Q.  Formally,  we  set 

Cm  =  {v  C  Q  |  v  =  uf=1  Si  where  Si, . . . ,  Si  C  Q  are  convex  and  closed}. 

Note  that  we  do  not  require  the  sets  Si, . . . ,  Si  to  be  distinct  so  that  C(ja  C  C(z),  for  any  k  <  i.  In  what  follows 
we  study  the  quotient  set  of  t- convex  and  closed  subsets  Cm  =  C(^)/~.  The  next  result  is  the  main  result  of 
this  section. 

Theorem  5.1  (Compactness  of  C^)  The  pair  (C^ ,  d@)  a  metric  space  and ,  with  the  topology  induced  by 
dg ,  the  set  C*^  is  compact. 


Proof.  It  is  easy  to  verify  that  dg  is  a  metric  on  C* .  Instead,  proving  the  compactness  of  Cm  requires  some 
attention.  We  aim  to  show  that  any  sequence  in  C*^  has  a  subsequence  that  converges  to  a  point  in  C*^ .  This 
fact’s  proof  is  articulated  in  several  steps  and  relies  upon  several  known  results: 


(i)  the  space  C,  endowed  with  the  Hausdorff  distance  dn  ■  C  X  C 
metric  space; 


^>o,  is  m  Theorem  0.8]  a  compact 


(ii)  if  a  sequence  of  closed  convex  subsets  of  Q  converges  in  the  Hausdorff  sense  to  a  set  K,  then  [2H  Propo¬ 
sition  1.6.8]  K  is  closed  and  convex;  and 

(iii)  for  any  two  convex  subsets  A,  B  of  Q  C  Rd,  it  is  known  [TS]  Eq.  (1)]  that 


de{A ,  B)  < 


/  2  hid 

(D' 

d- 1 

\  A 

\21/d  -  1 

l  2  , 

)  ) 

' )dH{A ,  B), 


where  D  =  max  {diam(A),  diam(B)}  and  where  Kd  is  the  volume  of  the  unit  sphere  in 


Now  let  {u*(n)}ngz> 0  be  any  sequence  in  Cm.  For  all  n  €  Z>o,  pick  any  representative  of  the  equivalence  class 
v*(n),  denote  it  by  v(n)  and  consider  the  sequence  {n(n)}raez>0-  Because  {v(n)}nez>0  is  a  sequence  in  and 
because  C  C,  it  follows  from  fact  (i)  above  that  {n(n)}nez>0  contains  a  subsequence  that  converges  in  the 
Hausdorff  sense  to  an  element  of  C.  In  other  words,  there  exist  {v(ni-)}k£i,>0  Q  {r(n)}„gz>0  and  v  £  C  such 
that 

lim  v(nk)  =  v, 

k — >oo 


where  (=  denotes  convergence  in  the  Hausdorff  sense.  We  claim  that  v  £  C(j)  so  that  the  set  Cm  is  compact.  To 
show  this  claim,  we  plan  to  exhibit  a  collection  of  convex  and  closed  subsets  of  Q ,  say  {Si, . . . ,  Sz},  such  that 


v  =  U  f=1Si. 


We  begin  to  prove  this  claim  as  follows.  By  definition  of  C^,  for  all  k  £  Z>o  there  exists  a  collection 
(Si(nfc), . . . ,  Sg(nk)},  of  convex  and  closed  subsets  of  Q  whose  union  equals  v(nk)-  Now,  we  consider  the 
sequence  {Siink)}k&>0-  Again,  since  (C,  dn)  is  a  compact  metric  space  we  have  that  there  exists  a  subsequence 


{fifci }fciez>0  Q  {nk}kez>0  such  that  linq-j^oo  Si(nkl)  ==  Si  for  some  Si.  Fact  (ii)  above  ensure  that  Si  is  a 
convex  closed  subset  of  Q.  Consider  now  the  sequence  {S2(nkl)}kl^>0.  By  reasoning  as  previously  we  know 

that  there  exists  a  subsequence  {nk2}k2 ez>0  Q  ink1}k1e'z>o  such  that  lim^^oo  S2(nk2)  (=  S2  where  S2  is  some 
convex  closed  set  of  Q.  Moreover,  it  is  clear  that  also  limfc2_>00  S\(nk2)  (=  Si.  By  iterating  this  procedure  we 
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conclude  that  there  exist  a  sequence  {ns}sez>0  and  a  collection  of  convex  and  closed  sets  {Si, . . . ,  Sf}  such  that 

Hindoo  Si(ns)  =  Si  for  all  i  £  { 1, . . .  ,£}. 

Next,  we  aim  to  show  that 


it  t 

lim  I  J  Si(ns)  =  I  J  lim  Si(ns )  =  I  J  St. 

s — >-oo  '  s — >-oo 

i— 1  2=  1  2=1 

For  simplicity,  let  us  denote  lim^oo  (Ji=1  Si(ns)  by  Soo-  For  p  £  Q1  note 


P  <£  Sa 


dist  ^p,  hm  (J  Si(ns)^  =  0 


lim  dist  I  p.  (J  Si{ns)  J  =  0, 


(15) 


where,  for  a  given  closed  set  X ,  dist(p,  X)  denotes  the  Euclidean  distance  between  p  and  X ,  namely,  minxex  ||p— 
cc|| .  Using  the  fact  that,  for  given  closed  sets  X  and  Y ,  dist(p,X  U  Y)  =  min  {dist(p,  X),  dist(p,  F)},  we  can 
write 


lim  dist 


P ,  U  Si(ns 

V  *=  l 


)  =0 


lim  min  {dist(p,  Si (ns)), ...,  dist(p,  S^(ns))}  =  0 

s — >oo 

3 j  £  (1, . . .  ,£}  s.t.  lim  dist(p,  S,  (ns))  =  0 

s — >oo 

3j  £  (1, s.t.  p  £  Sj 


t 

P  G  |J  St 
2=1 


The  above  chain  of  implications  proves  m ■  Now  observe  that,  from  the  uniqueness  of  the  limit  it  follows  that 
limfe^oo  v(nk)  =  lim  s — >oo  it  :1  Si(ns)  =  Ui=1  Si  =  v.  Since  Si  is  closed  and  convex  for  all  i  £  (1, . . .  ,£},  it 
follows  that  v  £  and,  in  turn,  that  C ^  endowed  with  the  Hausdorff  metric  is  a  metric  compact  space. 

To  establish  the  statement  of  the  Theorem  it  finally  remains  to  prove  that  either  lim/c_,00  dQ(v(rik),v)  =  0 
or  lirm^oo  de(i;(ns),  v)  =  0.  To  this  end,  observe  that,  given  X\,  X2,  Yi,  Y^  C  Q,  the  following  inclusion  holds 
(Xi  U  X2)  0  (Yi  U  Y2)  C  {X\  0  Fi)  U  (X2  0  Y2).  Hence,  we  compute 


v  ( ns )  ©  v 


which  implies 


(J  Si{ns) 


©  v  = 


e 

C  |J  (Si(ns)0Si), 
2=1 


de  ( v  {ns),v)  =  de 


t 

<J2de  (Si(n.),Si) 

i=l 


< 


< 


2  Kd  (  diam(<5) 

21/d  -  1  V  2 


{diam  (Sj(ns)) ,  diam  (Si) } 
2 

V  d—1  t 

■)  Y,dH  (Si(ns),Si), 

'  2=1 


d—1 

dH  (, Si(ns),Si ) 


where  the  second  and  the  third  inequalities  follow,  respectively,  from  fact  (iii)  above  and  from  the  upper  bounds 
max{diam  (Si  (ns)) ,  diam  (Si) }  <  diam(Q),  for  all  i  £  {1,  Since  lirm^oo  d#  (Si(ns),  S,;)  =  0  for  all 

i  £  we  conclude  that 


lim  de(v(ns),  v)  =  0. 


Now,  let  v*  denote  the  equivalence  class  for  which  v  is  a  representative.  Since  the  metric  cZq  is  insensitive  to  sets 
of  measure  zero,  it  follows  that  lim^oo  dQ(v* (na),v*)  =  0  and,  in  turn,  that  {v*(n)}nez>0  has  a  subsequence 
convergent  to  point  of  C*^.  This  concludes  the  proof  that  CX  is  a  compact  space.  I 
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(16) 


The  metric  de  naturally  extends  to  a  metric  over  the  space  (C*)N ,  and  hence  over  (Cfy)N ,  by  defining 

N 

de(u,v)  =  ^deiu^Vi), 

2=1 

for  any  u  =  {ut)f=1  and  v  =  (vi)fL i  in  (C*)N.  The  compactness  of  (Cfo)N  is  then  a  simple  consequence  of 
Theorem  15.11 

Corollary  5.2  (Compactness  of  (C^)N)  The  pair  ( (C  V ) ,V ,  d0 )  is  a  metric  space  and,  with  the  topology 
induced  by  dQ,  (C^)N  is  compact. 

6  On  the  continuity  of  some  geometric  maps  and  the  resulting  con¬ 
vergence  proofs 

In  this  section  we  prove  the  main  convergence  theorems  for  our  gossip  coverage  algorithms.  First,  however,  we 
need  to  establish  the  continuity  properties  of  certain  geometric  maps  and  of  the  proposed  algorithms  T  and  Ts . 
Before  proceeding,  we  discuss  two  significant  modeling  aspects.  First,  recall  that  Section[5]introduces  the  spaces 
C,  C(£),  CN ,  C*  and  the  corresponding  quotient  sets  C*,  C%y  (C*)N ,  (C^AN .  We  can  do  the  same  with  the 
partition  space  Vat.  Indeed  from  Definition^. l|[iii)|we  have  that  each  component  Vi  of  v  £  Vn  can  be  mapped  by 
the  canonical  projection  into  an  equivalence  class  v*  in  C*  \  {0};  in  turn,  any  v  £  Vn  can  be  naturally  associated 
to  the  IV-collection  of  equivalence  classes  v*  =  .  Accordingly,  we  denote  the  space  of  equivalence  classes 

of  N -partitions  by  Vn-  Recall  that  all  these  quotient  spaces  are  metric  spaces  with  respect  to  the  symmetric 
distance  dg. 

Second,  recall  the  following  maps:  the  centroid  map  Cd  :  {A  £  C  \  p(A)  >  0}  — >  Q  defined  in  equation  (©, 
the  1-center  function  Ti\  :  Q  x  C  — >  M>o  defined  in  equation  ©,  the  multicenter  function  7fCentroid  :  Vat  — >  K>o 
defined  in  equation  (Tl2l).  the  gossip  coverage  map  :  Vn  — >  Vat  with  i  y  j  defined  in  equation  (fill),  and, 
for  S  >  0  and  i  y  j,  the  modified  gossip  coverage  map  T A-  :  Vn  — ►  Vn  defined  in  Section  l3Al  We  claim  that 
all  these  maps  are  insensitive  to  sets  of  measure  zero.  To  substantiate  this  claim,  observe  that  the  integrals 
of  a  bounded  measurable  function  over  a  set  A  and  over  a  set  B  are  equal  if  de(A,  B)  =  0.  This  observation 
allows  us  to  redefine  the  centroid  map,  the  1-center  function  and  the  multicenter  function  as  Cd  :  C*  \  {0}  — >  Q, 
Hi  :  Q  x  C*  — >  R>o  and  Hcentroid  :  Vn  R>o>  respectively.  Regarding  the  modified  gossip  coverage  map  , 
we  reason  as  follows.  For  v*  £  V^,  let  v  £  Vn  and  v'  £  Vn  be  two  representatives  of  v*  and  let  v  and  v'  denote, 
respectively,  the  images  of  v  and  v'  under  the  map  T^,  that  is,  v  =  T^(v)  and  v'  =  T^{v').  Since  the  centroid 
map  and  the  definitions  of  the  points  pi  and  pj  in  equation  (1 1411  are  insensitive  to  sets  of  measure  zero,  it  follows 
that  d0  (v,vr)  =  0;  in  other  words  v  and  v'  belong  to  the  same  equivalence  class,  say  v*  £  From  these  facts 
we  can  redefine  the  modified  gossip  coverage  map  as  JA  :  Vn  — ►  An  analogous  argument  applies  to  the 

map  T.  This  concludes  the  justification  of  our  claim.  Finally,  note  that  the  Voronoi  map  V  :  QN  \  Sn  — >  Vn 
defined  in  equation  ©  can  be  composed  with  the  standard  quotient  projection  map  and  therefore  denoted  by 
V  :  Qn\Sn  —■ -  Viv-  summary,  the  centroid  map,  the  1-center  function,  the  multicenter  function,  the  modified 
gossip  coverage  map,  and  the  Voronoi  map  are  indeed  insensitive  to  sets  of  measure  zero. 

6.1  Continuity  of  various  geometric  maps 

We  start  by  recalling  that  the  compact  connected  set  Q  is  equipped  with  a  bounded  measurable  positive  function 
(j)  :  Q  — >  K>o-  We  define  the  diameter  of  Q  and  the  infinity  norm  of  </>  by  diam(Q)  =  max{||a;  —  y\\  \  x,y  £  Q} 
and  || 0|| oo  =  max^gQ  4>{x),  respectively.  The  following  lemma  states  some  important  properties  of  the  1-center 
cost  function. 

Lemma  6.1  (Continuity  properties  of  the  1-center  function)  Given  a  compact  convex  set  Q  C  Rd,  let 

(f>  :  Q  — >  R>o  be  bounded  and  measurable  and  let  f  :  R>o  — >  R>o  be  locally  Lipschitz,  increasing,  and  convex. 
Define  the  function  Hi  :  Q  x  C*  — >  R>o  as  in  equation  The  following  statements  hold: 

(i)  the  function  p  e- >  Hi(p;  A)  is  strictly  convex  in  p,  for  any  A  £  C*  \  {0}, 
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(ii)  the  function  p  i— >  7ii{p\  A)  is  globally  Lipschitz  in  p,  for  any  A  £  C* ,  and 
(Hi)  the  function  A  i— >  Tt\(p]  A)  is  globally  Lipschitz  in  A  with  respect  to  dQ,  for  any  p  £  Q. 
We  now  state  the  main  result  of  this  subsection. 


Theorem  6.2  (Continuity  properties  of  centroid,  multicenter  and  Voronoi  maps)  Given  a  compact  con¬ 
vex  set  Q  C  R'1,  let  <j>  :  Q  —>  R>o  be  bounded  and  measurable  and  let  f  :  R>o  — >  R>o  be  locally  Lipschitz, 
increasing,  and  convex.  With  respect  to  the  topology  induced  by  dQ,  the  following  maps  are  continuous: 


(i)  the  centroid  map  Cd  :  C*  \  {0}  — >  Q, 

(ii)  the  multicenter  function  7dCOntroid  :  — >  M>o, 

(Hi)  the  Voronoi  map  V  :  QN  \  Sn  — >  Vfj, 

(iv)  for  all  i,j  £  (1, . . . ,  TV},  i  b  j,  the  gossip  coverage  map  Tij  :  £  V%  |  Cd(^j)  b  Cd(nj)}  — >  Vfj,  and 

(v)  for  all  5  >  0,  i,j  £  (1, . . . ,  N},  i  b  j,  the  modified  gossip  coverage  map  J fj  :  Vfj  — >  Vfj. 


The  continuity  properties  (ii)  and  (iv)  (respectively,  (v) )  are  exactly  what  is  needed  to  apply  the  Krasovskii- 
LaSalle  invariance  principles  stated  in  Section  0]  to  the  gossip  coverage  algorithm  (respectively,  to  the  modified 
gossip  coverage  algorithm).  The  continuity  properties  (i)  and  (iii)  are  intermediate  results  of  independent 
interest. 

Proof.  [Proof  of  Lemma RT7T1  Let  Lf  be  the  Lipschitz  constant  of  /  :  [0,diam(Q)]  — >  R>o-  We  check  the  claims 
in  order,  beginning  with  statement  (i)  For  A  £  (0, 1),  A  £  C*  \  {0},  and  pi,p2  £  Q ,  we  compute 


Wi(Api  +  (1  -  X)p2-,A)  =  I  /(||Api  +  (1  —  A)p2  —  q\\)<P(q)dq 

J  A 

<  [  f(M\Pi  -  all  +  (1  -  A)||p2  -  q\\)(p(q)dq 

J  A 

<  JA  (A/(lbi  -  all)  +  (1  -  A)/(|b2  -  q\\j)<t>(q)dq 
=  XHi (pi ;  A)  +  (1  -  \)Hi(p2;A), 


(17) 

(18) 


where  inequality  ED  follows  from  the  triangle  inequality  and  from  /  being  increasing,  and  inequality  (11811 
follows  from  the  convexity  of  /.  This  inequality  proves  convexity.  Moreover,  since  the  first  inequality  is  strict 
outside  the  line  passing  through  p\  and  p2  and  since  A  has  non-empty  interior,  the  function  is  in  fact  strictly 
convex.  Note  that  statement  (i)  implies  that  p  >— ►  Hi(p;A)  is  locally  Lipschitz,  using 


Theorem  10.4].  The 


stronger  statement  (ii)  can  be  derived  as  follows.  For  pi,p2  £  Q ,  we  compute 


\7-li[pi,A) 


A)\  = 


< 

< 


I  f(\\pi-q\\)4>(q)dq-  f  f{\\P2-q\\)<t>(q)dq 

J  A  J  A 

I  [f(\\Pl  ~  9ll)  ^  f(\\P2  -  q\\)](f{q)dq 
J  A 

[  \f(\\Pi  ~  all)  -  f(\\P2  -  q\\)\<t>(q)dq 

J  A 

[  Lf\\p1-p2\\(j){q)dq  <  AHI^IUM^IIpi  -  p2||- 
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This  implies  the  Lipschitz  condition  in  statement  (ii)  Statement  (iii)  can  be  proved  as  follows.  Let  A,  A!  be 
two  elements  of  C* ,  note  A  =  {A  \  A')  U(A  fl  A')  and  compute 


\Hi{p\A)  -  7di(p;  A')\  = 


< 


[  f(\\p-q\\)<l>(q)dq 

-  1 

/(lb  -  q\\)<t>{q)dq 

JA\A' 

J 

A'\A 

1  f(\\p-q\\)<l>(q)dq 

+ 

[  f(\\p-q\\)<l>(Q)dq 

J  A\A' 

JA'\A 

I  f{\\P~Q\\)<l>{q)dq+  f 

/(lb-  q\\)<t>{<i)dq 

)  A\A' 

J  / 

V\A 

[  f(\\p-Q\\)<t>{q)dq 

)  AQ  A' 

<  max{/(||p- g||)  \  p,q  £  Aq  A'}  ||0||oo  ©  A') 

<  /(diam(Q))  ||<^||oo  ^e(A  A'), 


where  last  inequality  follows  from  /  being  increasing.  The  bound  implies  the  Lipschitz  condition.  I 

Before  proving  Theorem  16. 21  we  need  the  following  lemma  about  perturbations  of  convex  optimization  prob¬ 
lems. 

Lemma  6.3  Given  a  compact  convex  set  Q  C  and  a  metric  space  (X,d),  let  H  :  Q  x  X  — ►  R  have  the 
properties  that 

(i)  the  map  x  i— >  H(q,  x)  is  globally  Lipschitz  for  all  q  £  Q,  and 
(ii)  the  map  q  e- >  H(q,x)  is  continuous  and  strictly  convex. 

Then  the  map  q*  :  X  Q,  defined  by  q*(x)  =  argmin qGQ  H(q,x),  is  continuous. 


Proof.  Let  L h  be  the  Lipschitz  constant  of  x  >  H(q,x).  Thanks  to  the  Lipschitz  condition  of  the  function 
x  i— >  H(q,x),  for  all  a :,y  £  X,  the  point  of  minimum  q*(x )  takes  value  in  S'  =  {q  £  Q  \  H(q,y)  <  H(q*(y),y)  + 
2 L}jd{y,x)}.  Since  S'  is  a  sub-level  set  of  the  strictly  convex  function  q  i— ■>  H(q,y),  and  since  the  diameter  of 
a  sub- level  set  depends  continuously  on  the  level,  the  distance  ||  q*(y)  —  g*(x)l!  can  be  made  arbitrary  small  by 
reducing  d(y,x).  This  implies  the  claimed  continuity.  I 


We  are  now  ready  to  prove  the  main  result.  Proof.  [Proof  of  Theorem  16. 2j  We  prove  the  theorem  claims  in 

the  order  in  which  they  are  presented.  Claim  (i)  follows  combining  Lemma  16.31  and  Lemma  16.11  Since  the 
multicenter  function  7iCentroid  is  a  sum  of  suitable  1-center  functions  TL\,  the  claim  (ii)  is  also  immediate. 

Regarding  claim  (iii)  we  discuss  in  detail  the  two  dimensional  case.  We  first  prove  the  continuity  when 
N  =  2:  let  pi  and  p2  be  points  in  Q.  Let  21  =  ||pi  —  P2\\-  Since  pi  ^  p2,  l  >  0.  Up  to  isometries,  we  can 
assume  that,  in  the  Euclidean  plane  (x:y),  p\  =  (—1, 0)  and  p2  =  (^,0).  Let  d\  and  e?2  be  the  distances  from 
the  origin  of  points  pi  and  p2,  respectively.  It  is  clear  that  the  two  Voronoi  regions  of  pi  and  p2  are  separated 
by  the  locus  of  points  {x  £  Q  \  ||x  —  pi||  =  ||cc  —  P2II},  that  is  the  vertical  axis.  Now,  we  assume  that  the 
positions  of  p\  and  P2  are  perturbed  by  a  quantity  less  than  or  equal  to  e,  with  0  <  e  <  l.  By  effect  of  the 
perturbation,  the  axis  separating  the  two  Voronoi  regions  is  perturbed,  but  it  is  contained  in  the  locus  of  points 
Y12 (e)  =  {x  £  Q  |  ||a;  -piH  —  || z  —  P2II  <  2e}.  By  definition,  this  is  the  set  comprised  between  the  two  branches 
of  the  hyperbola  whose  equation  is  t-  -  rrin  =  1-  By  elementary  geometric  considerations,  the  area  of  this 
region  can  be  upper  bounded  by 


li(Y12(e))  <2e2  diam(<5)  +  4diam(Q)2 — ' ^  <  4  diam(<5)fl  + 

1  1  ' 


This  bound  implies  the  continuity.  The  case  in  which  N  >  2  follows  because,  moving  all  points  by  at  most  e, 
the  change  in  all  the  regions  is  upper  bounded  by  Ui<i  j<N  which  vanishes  as  e  — >  0+. 
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The  last  remaining  step  is  to  prove  claims  (iv)  and  (v)  We  focus  on  claim  (v)  and  show  claim  (iv)  as  a 


byproduct.  Let  v  G  Vn  and  let  v  =  T^{v).  According  to  the  definitions  in  Subsection  13.41  v  is  characterized  by 
the  sets  Ri  and  Rj .  Recall  that  these  two  sets  depend  on  the  sets  Ri  Rj ,  on  the  scalar  (v)  and  on  the  points 
Pi  and  pj.  One  can  see  that  Pij(v)  is  a  continuous  function  of  its  arguments  and  Vj.  Hence,  it  suffices  to 
show  that  also  Ri,  Rj  and  pi,  pj  depend  continuously  on  Vi  and  Vj .  To  do  this,  introduce  v'  G  Vn  and  compute 
v'  =  Tfj( v').  Assume  Cd(wj)  ^  Cd(i>j)  and  Cd(w')  ^  Cd(u').  Analogously  to  how  we  defined  Ri  and  Rj,  we  now 
define  the  regions  R!i  =  v[  fl  -ffbisector(Cd(u'-),  Cd('u'))  and  R'-  =  u'  n  Hbisector(Cd(u'),  Cd(w')).  We  aim  to  upper 
bound  the  composite  distance  +  d0 (Rj ■,  Rj)-  Observe  that  this  composite  distance  depends  on  the 

difference  of  the  two  argument  regions  Vi,v\  and  Vj ,  v\  both  directly  and  indirectly  via  the  induced  difference 
between  the  centroids.  Recalling  the  proof  of  claim  (iii)  let  e  be  an  upper  bound  on  the  displacement  between 
the  two  centroids.  Then  the  region  Y(e)  =  {x  G  Q  | ||rc  —  Cd(ui)||  —  ||x  —  Cd(wj)||  <  2e}  needs  to  be  included 
in  the  upper  bound  on  the  composite  distance.  Combining  these  considerations  we  obtain 


de{Ri,R'i)  +  de(Rj,R'j)  <  (de(vi,Vi)  +de(vj,v'j)) 

+  /x(F(max{ 


Cd(«<)  -  CdKjll,  II  Cdfe)  -  Cd(w'.)ll}))-  (19) 


Clearly,  if  d0(uj,i;')  - 
turn,  also  d0(i?i,I?') 

the  regions  Vi  and  Vj .  This  is  enough  to  prove  statement  (iv)  provided  the  two  regions  Vi  and  Vj  have  distinct 


0  and  de(vj,Vj)  — >  0,  then  ||  Cd(uj)  —  Cd(t>')||  — >  0  and  ||  Cd(uj)  —  Cd(i;')||  — >  0  and,  in 
dQ(Rj,  Rj)  — *  0.  Hence,  we  can  argue  that  the  sets  Ri  and  Rj  depend  continuously  on 


centroids.  Moreover  a  direct  consequence  of  this  fact  is  that  also  the  points  pi  and  pj  depend  continuously  on 
Vi  and  Vj.  Finally,  observe  that  in  the  limit  case  Cd(ui)  =  Cd(u,-)  the  continuity  of  Ri  and  Rj  is  captured  by 
the  fact  that  Pij(v)  is  a  continuous  function  of  Cd(vj), Cd(i>j)  and  that  Pij(v)  =  0  if  Cd(vi)  =  Cd(uj).  The 
continuity  of  Ri,  Rj,  Pij(v),  Pi  and  Pj  imply  also  the  continuity  of  Ri  and  Rj  and,  in  turn,  of  TA.  I 


6.2  Convergence  proofs 

In  view  of  the  identification  between  ^partitions  and  their  equivalence  classes  introduced  at  the  beginning  of 
this  section,  we  are  now  ready  to  complete  the  proof  of  the  convergence  results  presented  in  Section  13.31 

We  start  by  clarifying  the  precise  meaning  of  convergence  in  Theorems  13.71  and  13. 81  Specifically,  we  say  that 
a  sequence  of  partitions  {v(t)}t^z>0  C  Vat  converges  to  a  set  of  partitions  X  C  Vn  if  the  symmetric  distance 
from  {v(t)}t<zZ> o  to  X  converges  to  zero,  that  is, 

lim  inf{d0(?;(t), x)  \  x  G  X}  =  0. 


Proof.  [Proof  of  Theorem  13. 8|  We  prove  the  deterministic  statement  (i).  We  start  by  observing  that,  through 

the  canonical  projection,  the  evolution  {v(t)}tG%>0  C  Vn  of  Ts  :  Vn  — ►  Vn  can  be  mapped  into  the  evolution 
{i>*(t)}teZ>0  c  of  Ts  :  We  aim  to  apply  Theorem  14. 31  to  the  dynamical  system  Ts 

and  its  evolution  {v* (t)}tez>0  C  V^.  In  what  follows,  our  goal  is  to  verify  whether  Assumptions  (i) 
and  (iv)  of  Theorem  14.31  are  satisfied. 


v; 


V 


Mil  (iii) 


Since  {v(t)}t^z>0  is  non-vanishing  and  finitely  convex  by  assumption,  it  follows  that  there  exists  <Gtf  such 
that  the  H-limit  set  of  {v(t)}tez>0  is  contained  in  fl  Vtv,  that  is,  fi(u(t))  C  C(^  fl  Vat.  This  implies  also 

that  fl(v*(t))  C  CX  fl  Vff.  As  stated  in  Theorem  15.11  CX  is  compact  in  the  topology  induced  by  the  metric 


Pq.  Hence,  even  though  CX  is  not  strongly  positive  invariant  for  T°,  the  weaker  version  of  Assumption  (i)  of 
Theorem  m  as  given  in  Remark R~TH  holds  true  for  the  sequence  {v* (t)}teZ>0  G  V^.  Now,  as  one  can  deduce 
from  Theorem  I6.^1ii)|  and  Lemma  13.41  the  function  7Ycentroid  satisfies  the  Assumption  (ii)  of  Theorem  14.31 
thus  playing  the  role  of  a  Lyapunov  function  for  the  dynamical  system  Ts . 


Moreover,  from  Theorem  I6.^(v)[ 
note  that  the  system  evolves  through  maps  that  are  continuous  in  Vfj  with  respect  to  the  metric  thus 
the  Assumption  (iii)  of  Theorem  14.31  is  satisfied.  Finally  observe  that  the  Assumption  (iv)  of  Theorem  14.31 


corresponds  to  the  assumption  of  uniform  persistency  in  Theoreml3.8l  Therefore,  we  conclude  that  the  evolution 
{v*(t)}tez> o  converges  to  the  intersection  of  the  fixed  points  of  the  maps  T^,  for  all  i,j  G  {1, . . . ,  A^},  j  ^  i. 
According  to  Lemma  |3. 31  this  intersection  coincides  with  the  set  of  mixed  centroidal  Voronoi  partitions  up  to 
sets  of  measure  zero. 
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The  proof  of  the  stochastic  statement  (ii)  follows  the  same  lines,  applying  Theorem  14.51  instead  of  Theo¬ 
rem  14.31  ■ 


Proof.  [Proof  of  Theorem  13. 7|  The  proof  of  this  result  follows  the  lines  of  the  proof  above  with  two  distinctions. 
The  distinctions  come  from  the  assumption  that  the  evolution  v,  in  addition  to  being  non- vanishing  and  finitely 
convex,  is  also  distinct  centroidal.  The  first  distinction  is  as  follows.  In  order  to  apply  the  Krasovskii-LaSalle 
invariance  principle  we  require  the  continuity  property  stated  in  Theorem  IG.^liv)  Additionally,  we  note  that 
the  space  of  finitely  convex  and  distinct  centroidal  partitions 


{v  G  Vn  |  ||  Cd(ui)  —  Cd(i>j)||  >  e  for  all  i  ^  j} 

is  a  closed  and  hence  compact  subspace  of  Vat.  The  second  distinction  is  as  follows.  Since  we  rule  out  the  case  of 
coincident  centroids,  we  can  infer  convergence  to  centroidal  Voronoi  partitions  instead  of  convergence  to  mixed 
centroidal  Voronoi  partitions;  see  Lemma  13. t|fii)|  ■ 


7  Conclusions 

In  summary,  we  have  introduced  novel  coverage,  deployment  and  partitioning  algorithms  for  robotic  networks 
with  minimal  communication  requirements.  To  analyze  our  proposed  algorithms,  we  have  developed  and  char¬ 
acterized  (1)  intuitive  versions  of  the  Krasovskii-LaSalle  Invariance  Principle  for  deterministic  and  stochastic 
switching  systems,  (2)  relevant  topological  properties  of  the  space  of  partitions,  and  (3)  useful  continuity  prop¬ 
erties  of  a  number  of  geometric  and  multicenter  functions. 

We  believe  there  remain  interesting  open  issues  in  the  study  of  gossiping  robots  and  of  dynamical  systems  on 
the  space  of  partitions.  We  are  keen  on  extending  these  ideas  to  non-convex  complex  environments  and  discrete 
environments  such  as  graphs.  Following  Remark  13.21  we  plan  to  study  gossip  coverage  algorithms  for  more 
general  multicenter  functions,  including  nonsmooth,  anisotropic  and  inhomogeneous  functions.  Additionally,  we 
plan  to  investigate  gossip  coverage  algorithms  capable  of  adapting  to  time- varying  scenarios  such  as  problems  in 
which  robotic  agents  arrive  to  and  depart  from  the  network.  Finally,  inspired  by  stigmergy  in  territorial  animals, 
we  plan  to  design  communication  protocols  for  multiagent  systems  based  on  the  ability  of  leaving  messages  in 
the  environment. 


A  A  robotic  network  implementation  of  gossip  coverage  algorithms 

In  Subsection  cm  we  discussed  coverage  control  algorithm  for  groups  of  robots  with  synchronized  and  reliable 
communication  along  all  edges  of  a  Delaunay  graph;  then  in  Subsection  13.21  we  introduced  our  gossip  coverage 
algorithms.  Here  we  discuss  in  full  detail  one  possible  way  of  implementing  the  partition-based  gossip  coverage 
algorithm  in  a  robotic  network  with  weak  communication  requirements. 

We  consider  a  group  of  agents  all  having  the  following  capabilities:  (Cl)  each  agent  *  G  {1, . . . ,  N}  knows  its 
position  and  moves  at  positive  speed  iq  to  any  position  in  the  compact  convex  environment  Qcl2;  (C2)  each 
agent  may  store  an  arbitrary  number  of  locations  in  Q  and  has  a  clock  that  is  not  necessarily  synchronized  with 
other  agents’  clocks  (specifically,  we  assume  same  clock  skew,  but  different  clock  offsets  among  the  agents);  and 
(C3)  if  any  two  agents  are  within  distance  rComm  of  each  other  for  some  positive  duration  of  time,  then  they 
exchange  information  at  the  sample  times  of  a  Poisson  process  with  intensity  ACOmm- 

The  random  destination  &  wait  algorithm  is  described  as  follows.  Given  a  parameter  e  <  rComm/4,  each  agent 
i  G  (1, . . . ,  N}  maintains  in  memory  a  dominance  region  Vi  and  determines  its  motion  by  repeatedly  performing 
the  following  three  actions  over  periods  of  time  that  we  label  epochs: 

1:  instantaneously  agent  i  selects  uniformly  randomly  a  destination  point  qi  in  the  set  {q  G  R2  |  dist  ( q,dvi  \ 

dQ)  <  e}; 

2:  agent  i  moves  in  such  a  way  as  to  reach  point  qi  in  time  precisely  equal  to  d  =  diam(Q)/  min{iii, . . . ,  ujv}; 
and 

3:  agent  i  waits  at  point  qi  for  a  duration  that  equals  either  d  with  probability  1/2  or  2d  with  probability  1/2. 
We  have  assumed  that  agents  may  reach  locations  at  a  distance  up  to  e  away  from  Q\  this  assumption  can  be 
removed  at  the  cost  of  additional  notation. 
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The  random  destination  &  wait  algorithm  is  meant  to  be  implemented  concurrently  with  the  modified 
gossip  coverage  algorithm  -  where  the  parameter  <5  is  selected  to  satisfy  <5  <  rCOmm/4.  The  two  algorithms 
jointly  determine  the  evolution  of  the  agents  positions  and  dominance  regions  as  follows.  If  any  two  agents 
i  and  j  are  within  communication  range  rCOmm  at  any  instant  of  time  and  for  some  positive  duration,  then, 
at  each  sample  time  of  the  corresponding  communication  Poisson  process,  the  two  agents  exchange  sufficient 
information  to  update  their  respective  regions  i’i  and  Vj  via  the  modified  gossip  coverage  map  T^. 

Proposition  A.l  (Random  destination  &:  wait  ensures  persistent  random  gossip)  Consider  a  group 
of  agents  with  capacities  (Cl),  (C2)  and  (C3)  and  parameters  Ui,  rcomm  and  ACOmm-  Assume  the  agents  im¬ 
plement  the  random  destination  &  wait  algorithm  and  the  modified  gossip  coverage  algorithm  with  parameter 
e  <  rc omm/4  and  S  <  rCOmm/4.  The  following  statements  hold: 

(i)  the  sequence  of  applications  of  the  modified  gossip  coverage  map  is  a  randomly  persistent  stochastic  process; 
and 

(ii)  the  resulting  evolution  v  :  R>o  — >  Vn ,  conditioned  upon  being  non-vanishing  and  finitely  convex,  converges 
almost  surely  to  the  set  of  mixed  centroidal  Voronoi  partitions. 


Proof.  Assume  that  at  time  s  £  R>o,  the  regions  vt  and  Vj  are  at  distance  6  or  less.  Pick  any  two  points 
Pi  £  dvi  \  dQ  and  pj  £  dvj  \  dQ  such  that  \\pi  —  Pj\ \  <26  and  define  two  balls  of  radius  e  centered  at  pi  and  pj 
by  Bi  =  {q  £  M2  |  \\q  —  Pi\\  <  e}  and  Bj  =  {q  £  !2  |  ||q  —  pj\\  <  e}.  By  construction,  each  point  in  Bi  is  at  most 
at  distance  2e  +  26  <  rcomm  from  each  point  in  Bj.  We  claim  that,  independently  of  the  agent  positions  and 
the  environment  partitions  at  time  s  and  of  their  evolution  throughout  the  past  time  interval  [0,  s],  there  exists 
a  duration  A  £  R>o  such  that  the  meeting  event  “agent  i  is  in  Bi  and  agent  j  is  in  Bj  during  the  time  interval 
[s  T  A,  s  +  A  +  g?]”  has  probability  that  is  lower  bounded  by  a  positive  constant.  If  these  meeting  events  have 
positive  probability,  then  communication  events  have  positive  probability  and  therefore  the  proposition  follows. 

We  prove  the  claim  as  follows.  At  time  s,  we  say  that  an  agent  is  in  OFF  state  if  it  is  at  the  initial  time  of  an 
epoch  or  it  is  moving  to  its  next  destination  point.  We  say  it  is  in  ONi  state  if  it  has  just  reached  a  destination 
point  or  it  has  been  waiting  at  the  destination  point  for  a  duration  of  time  less  than  d.  We  say  it  is  in  ON2 
state  if  it  has  been  waiting  at  the  destination  point  for  a  duration  of  time  greater  than  or  equal  to  d  and  less 
than  2d.  Note  that  each  agent  remains  in  each  one  of  the  three  states  {OFF,  ONi,  ON2}  for  a  duration  of  time 
d.  The  transitions  between  states  are  regulated  by  a  Markov  chain  with  three  states  {OFF,  ONi,  ON2}  and  with 
(column-stochastic)  transition  matrix 


1/2  1' 
0  0 
1/2  0 


For  now,  suppose  that  at  time  s,  agent  i  is  at  the  initial  time  of  an  interval  of  duration  d  in  either  state  OFF,  or 
state  ONi  or  state  ON2.  From  the  matrices 


'1/2 

1/4 

1/2 

'1/2 

3/8 

1/4" 

A3  = 

1/2 

1/2 

0 

,  and  A5  = 

1/4 

1/2 

1/2 

_  0 

1/4 

1/2 

1/4 

1/8 

1/4 

we  establish  the  following  three  facts: 

(i)  if  agent  i  is  at  state  OFF  at  time  s,  then  at  time  s  +  3d  it  is  at  state  OFF  with  probability  1/2  and  at  time 
s  +  5d  it  is  at  state  ON2  with  probability  1  / 4; 

(ii)  if  agent  i  is  at  state  ONi  at  time  s,  then  at  time  s  +  3d  it  is  at  state  OFF  with  probability  1/4  and  at  time 
s  +  5d  it  is  at  state  0N2  with  probability  1/8;  and 

(iii)  if  agent  i  is  at  state  ON2  at  time  s,  then  at  time  s  +  3d  it  is  at  state  OFF  with  probability  1/2  and  at  time 
s  +  5d  it  is  at  state  ON2  with  probability  1  /4. 


We  now  return  to  the  two  agents  i  and  j.  According  to  the  robot  capability  (C2),  at  time  s  the  clocks  of 
agent  i  and  agent  j ,  denoted  by  ti(s)  and  tj(s)  respectively,  are  distinct  and  can  be  written  as  an  integer 
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number  of  durations  d  plus  a  quotient;  specifically  U(s)  =  kid  +  rji ,  and  tj(s)  =  kjd  +  rjj,  with  ki,kj  G  Z >o , 
771 , 772  G  [0,g?[,  and  possibly  ki  7^  kj  and  77,  7^  7^-.  Without  loss  of  generality,  assume  rji  <  rjj.  Facts  (i),  (ii), 
and  (iii)  and  the  assumption  that  clock  skew  is  equal  among  agents,  jointly  imply  that,  with  probability  greater 
or  equal  to  (1/8)2  =  1/64,  agents  i  and  j  during  the  interval  [s  +  Ad  —  rji,  s  +  5d  —  rji]  are  waiting  at  some 
destination  point  chosen,  respectively,  at  time  s  +  3d  —  rji  and  s  +  3d  —  rjj .  If  agent  i  is  in  Bi  and  agent  j  is 
inside  Bj  during  the  interval  [s  +  Ad  —  rji,  s  +  5d  —  rji],  then  they  communicate  at  least  once  with  probability 
Pconm  :=  1  —  e-dAcomm  >  0.  Let  us  now  find  a  lower  bound  on  the  probability  of  agents  i  being  in  Bi  at  time 
s  +  Ad  —  rji]  clearly  this  happens  if  the  following  two  facts  occur: 

•  during  the  interval  [s,  s  +  3d  —  rji]  agent  i  does  not  communicate  with  any  other  agent  and  hence  its 
dominance  region  remains  unchanged;  and 

•  at  time  instant  s  +  3d  —  rji  agent  i  is  at  the  beginning  of  a  new  epoch  (state  OFF )  and  selects  a  destination 
point  in  Bi. 

Since  during  the  interval  [s  —  rji,  s  +  3d  —  rji]  agent  i  spends  a  time  duration  of  at  most  2d  at  distance  less 
than  S  from  any  other  agents,  then  the  probability  that  agent  i  does  not  communicate  in  that  interval  of 
time  is  lower  bounded  by  (1  —  pComm)2  =  e_2rfAcomm.  Moreover,  provided  that  at  time  instant  s  +  3d  —  rji 
agent  i  is  at  the  beginning  of  a  new  epoch  (state  OFF)  and  that  the  vt  is  unchanged  until  time  instant  s, 
we  have  that  the  destination  point  for  agent  i  is  selected  uniformly  randomly  so  that  it  belongs  to  Bi  with 
probability  area(-Bi)/ area  ({g  G  K2  |  dist  (q,  dvi  \  dQ)  <  e})  >  pare a  :=  7re2/(  area(Q)  +  e  perimeter(Q)) . 
Similar  considerations  hold  also  for  agent  j . 

In  summary,  if  agents  i  and  j  are  neighbors  at  time  s,  then  they  communicate  at  least  once  during  the  interval 
[s  +  Ad  —  rji,  s  +  5d  —  rji]  with  probability  lower  bounded  by  (l/64)p2rea(l  —  pComm)2Pcomm  >  0.  If  agents  i  and 
j  are  at  distance  greater  than  S  at  time  s  (and  therefore  they  are  not  neighbors),  then  the  modified  gossip 
coverage  map  is  the  identity  map  so  that  their  communication  is  immaterial.  I 


B  A  counterexample  showing  the  necessity  of  uniformly  persistent 
switches 


Theorem  14.31  contains  a  persistent  switching  conditions,  that  is,  it  requires  the  existence  of  D  G  N  such  that 
every  map  T),  i  G  {1, . . . ,  m},  is  applied  at  least  once  within  every  interval  [n,  n  +  D[,  n  £  Z>o-  This  appendix 
contains  an  example  proving  the  necessity  of  this  condition. 

Consider  the  plane  in  polar  coordinates  X  =  M>o  x  [0,  27r[U{0,  0}.  Define  the  standard  metric 
R>o  as  follows:  let  (pi,0i),  (P2A)  be  any  pair  of  elements  of  X  and 

d((pi,  0i),  {P2,02))  =  y/ (pi  cos  61  -  P2  cos  62)2  +  (pi  sindi  -  p2  sinfc)2. 

Consider  now  the  continuous  maps  T,  :  X  — ^  X,  i  G  {1,2},  defined  by  respectively 


Ti(p,0) 


T2(p,0) 


(P2,0), 

(^i,  (0  +  p-l) 

((1  -  sin  9)p,  6) , 
(P,0) 


if  0  <  p  <  1, 
mod  27t^  ,  if  p  >  1, 

if  0  <  9  <  7T, 
if  7 t  <  9  <  2tt. 


Define  T  :  X  =4  X  by  T(p,9)  =  {T\  ( p,9 )  ,T2  (p,9)}  and  the  function  U  :  X  — >  R>0  by  U  ( p,9 )  =  p.  Observe 
that  U  is  continuous  and  non-increasing  along  T.  Assume  now  that  there  exists  D  G  N  such  that,  for  any 
n  G  Z>o,  there  exist  ni  and  n2  within  the  interval  ]n,n  +  D]  such  that  xni+i  =  T\(xni)  and  xn2+i  =  T2(xn2). 
Then,  by  Theorem  14. 31  the  w-limit  set  of  each  evolution  of  T  is  a  subset  of 


{(p,  9)  G  X  |  p  =  1, 7T  <  9  <  2tt}  U  {0,  0}. 


(20) 


Next,  we  relax  the  condition  that  the  map  T2  is  applied  at  least  once  inside  each  interval  of  arbitrary  ampli¬ 
tude  D  and  we  show  that  there  exists  one  sequence  that  does  not  converge  to  the  w-limit  set  in  equation  (l20l). 
To  this  aim,  assume  the  sequence  {(p(n),  9(n))}nez>0  satisfies 
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(i)  p(0)  >  1; 

(ii)  (p(l),0(l))=T1(p(O),0(O))  and 


(iii)  {p(n+ 1),  9(n+ 1))  =  T2  ( p{n),9{n ))  if  and  only  if  n  <  9{n)  <  2n  and  (p(n),  9(n))  =  T\  ( p[n  —  1),  9{n  —  1)). 

Note  that  if  tv  <  9(n)  <  2n,  then  T2  ( p(n),9(n ))  =  ( p(n),9(n )).  Therefore,  the  evolution  {(p(n),  9(n))}  equals 
|(p(n),0(n))|  where  (p(O),0(O))  =  (p(0),9(0))  and  (p(n),  9(n))  =  T"  ^p(O),  9(0)^j.  Regarding  this  new  se¬ 
quence,  observe  that 


1  <  p(i)  <  2  and  p[i  +  1)  <  p{i),  for  all  i  >  1, 

0  <  9{i  +  1)  —  9(i)  <  7 r,  for  all  i  >  1,  and  lim  (d(i  +  1)  —  9(i)  \  =  0 

i — >oo  V  / 

r  r  r 

lim  y  ( 9(i  +  1)  —  0(*)')  =  lim  W  ( p{i )  —  1)  =  lim  W  (  ~ - f-  *  —  1  )  =  oo, 

r — »oo  '  J  \  /  r — ►oo  '  ^  i — >oo  '  ^  1 


i—  1 


i—  1 


feUD-i 


(21) 

(22) 

(23) 


where  the  equality  p(«)  —  1  =  +  i  —  1)  1  can  be  proved  by  induction  over  i.  Properties  (12T1).  (|22l). 

and  (l23|  ensure  that  there  exists  a  sequence  {rih  \  h  G  Z>0}  such  that  ( p{rih ),  9{rih ))  =  T2  (p(rih  —  1),  9{rih  —  1))) 
for  all  h  G  Z>o,  and  ( p(t),0(t ))  =  Ti  (p(t  —  1),  9(t  —  1)))  if  t  ^  {n^  |  h  G  Z>o}.  Moreover,  we  have  that 
lim^_>00(ri7l+i—  rih)  =  oo.  In  other  words,  both  the  maps  T\  and  T2  are  applied  infinitely  often  along  the  evolution 
described  by  {(p(n),  0(n))},  but  there  does  not  exists  DgN  such  that  T2  is  applied  at  least  once  within  each 
interval  [n,  n+D],  n  G  Z>0.  Observe  that,  in  this  case,  ( p(n),9(n ))  converges  to  the  set  {(1, 9)  \  9  £  [0,  2tt]}  C  X. 
This  set  is  different  from  the  w-limit  set  in  equation  (121)1) . 


C  Discontinuity  of  the  multicenter  function  in  the  Hausdorff  metric 


To  see  that  7iCentroid  :  Vfy  — >  R>o  is  not  Hausdorff-continuous,  consider  the  sequence  of  2-partitions  {v(t)}tez>0 
of  the  interval  [—1,1]  CR  defined  by 


Mt) 


1 

2^+1 


2t_1  —  1 

u  u 

h=—  (2t  — 1  —  1) 


h 


1  h 

2t+1  ’  2t_1 


2t+i 


1 

u 

1  2*+1’1_ 

and  by  v2 (i)  =  [ — 1, 1]  \  'i’i(f).  Note  that  both  sequences  {v\{t)}t^z>0  and  {v2(t)}t£z>0  converge  to  [—1,1],  and 
that  Cd(tq(i))  =  Cd (v2(t))  =  Cd([—  1, 1])  =  0,  for  all  t  >  0.  Hence,  for  <j)(s)  =  1  and  /(s)  =  s,  we  compute 

7ii(0,  vi(t))  =  I  \x\dx  =  2  I  \x\dx  =  2—  =  l, 

J v\ (t)  J v\ (t)  n[o,i]  2 

and  consequently  Wcentroid(^(£))  =  2,  while  Wcentroid(limi-KX)  v(t))  =  27^i(0,  [— 1, 1])  =  4.  This  shows  that 

lim^—^oo  7^centroid(^(^))  7^  ^centroid (HlUt— >00  ^(^))‘ 
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